♦ Human Brain Mapping 33:2005-2034 (2012) ♦ 



Probabilistic Tractography Recovers a 
Rostrocaudal Trajectory of Connectivity Variability 
in the Human Insular Cortex 

Leonardo Cerliani, 1 ' 2 ' 3 * Raj at M. Thomas, 4 Saad Jbabdi, 5 

12 6 I 2 12 7 

Jeroen C.W. Siero, ' ' Luca Nanetti, ' Alessandro Crippa, ' 
Valeria Gazzola, ' ' Helen D'Arceuil, and Christian Keysers ' 

1 BCN Neurolmaging Center, University of Groningen, A. Deusinglaan, 2-9713AW Groningen, 

The Netherlands 

2 Department of Neuroscience, University Medical Center Groningen, A. Deusinglaan, 
2-9713AW Groningen, The Netherlands 
3 Social Brain Laboratory, Netherlands Institute for Neuroscience, Royal Netherlands Academy of Arts 

and Sciences, Amsterdam, The Netherlands 
i Kapteyn Astronomical Institute, Faculty of Mathematics and Natural Sciences, 
University of Groningen, 9700AV Groningen, The Netherlands 
5 Oxford Centre for Functional Magnetic Resonance Imaging of the Brain (FMRIB), 
Department of Clinical Neurology, University of Oxford, United Kingdom 
h Department of Radiology, University Medical Center Utrecht, Utrecht, The Netherlands 
7 Scientific Visualization and Computer Graphics Group, Institute for Mathematics and Computing 

Science, University of Groningen, PO Box 800, 9700AV Groningen, The Netherlands 
^Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts 
General Hospital, Harvard Medical School, Charlestown, MA 02129 



♦ ♦ 

Abstract: The insular cortex of macaques has a wide spectrum of anatomical connections whose distri- 
bution is related to its heterogeneous cytoarchitecture. Although there is evidence of a similar 
cytoarchitectural arrangement in humans, the anatomical connectivity of the insula in the human brain 
has not yet been investigated in vivo. In the present work, we used in vivo probabilistic white-matter 
tractography and Laplacian eigenmaps (LE) to study the variation of connectivity patterns across insu- 
lar territories in humans. In each subject and hemisphere, we recovered a rostrocaudal trajectory of 
connectivity variation ranging from the anterior dorsal and ventral insula to the dorsal caudal part of 
the long insular gyri. LE suggested that regional transitions among tractography patterns in the insula 
occur more gradually than in other brain regions. In particular, the change in tractography patterns 
was more gradual in the insula than in the medial premotor region, where a sharp transition between 
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different tractography patterns was found. The recovered trajectory of connectivity variation in the 
insula suggests a relation between connectivity and cytoarchitecture in humans resembling that previ- 
ously found in macaques: tractography seeds from the anterior insula were mainly found in limbic 
and paralimbic regions and in anterior parts of the inferior frontal gyrus, while seeds from caudal in- 
sular territories mostly reached parietal and posterior temporal cortices. Regions in the putative dys- 
granular insula displayed more heterogeneous connectivity patterns, with regional differences related 
to the proximity with either putative granular or agranular regions. Hum Brain Mapp 33:2005-2034, 

2012. © 2011 Wiley Periodicals, Inc. 

Keywords: DTI; insula; laplacian eigenmaps; diffusion-weighted imaging; connectivity-based 
parcellation 
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INTRODUCTION 

Two major features characterize the anatomical organiza- 
tion of the insular cortex: the wide spectrum of its anatomi- 
cal connections and its heterogeneous cytoarchitecture 
[Mesulam and Mufson, 1985]. In macaques, the insula proj- 
ects to primary and secondary sensory cortices [Friedman 
et al., 1986; Mesulam and Mufson, 1982b; Mufson and 
Mesulam, 1982], to premotor regions [Luppino et al., 1993; 
Matelli et al., 1986], to the posterior parietal lobule [Ander- 
sen et al., 1990; Cavada and Goldman-Rakic, 1989; Rozzi 
et al., 2006] as well as to higher-order association areas in 
the prefrontal lobe [Petrides and Pandya, 1999]. The insula, 
especially its anterior portion, also has extensive connec- 
tions with the amygdaloid complex [Amaral and Price, 
1984; Mufson et al, 1981; Stefanacci and Amaral, 2002] as 
well as with other limbic and paralimbic regions in the 
frontal and temporal lobes [Cavada et al., 2000; Kondo 
et al, 2005; Mesulam and Mufson, 1985; Moran et al., 
1987]. Insular projections were found to both limbic and 
motor compartments of the striatum [Chikama et al., 1997; 
Fudge et al., 2005] to autonomic nuclei in the brainstem 
and several thalamic nuclei (see Augustine [1985, 1996] 
and Mesulam and Mufson [1985] for extensive reviews). 

The distribution of the insular cytoarchitectonics is 
organized radially around the limen insulae [Mesulam 
and Mufson, 1982a; Nieuwenhuys et al., 2008]: the anterior 
ventral part of the insula, coextensive with the piriform 
cortex and the claustrum, has an agranular structure, 
while the posterior dorsal part has an isocortical homotyp- 
ical structure. The transitional dysgranular cortex in the 
anterior-dorsal, middle, and posterior-ventral insular ter- 
ritories shows a continuous pattern of progressive lamina- 
tion with an incipient presence of granule cells in layers 
II and IV in the anteroventral-posterodorsal direction 
[Carmichael et al, 1994; Mesulam and Mufson, 1982a]. A 
similar architectural organization was also found in the 
human insula [Bonthius et al., 2005; Mesulam and Mufson, 
1982a; Rose, 1928; Shaw et al, 2008; Zilles, 2004] including 
recent observer-independent probabilistic cytoarchitectonic 
maps of the posterior insula [Kurth et al., 2010a]. 

In macaques' insula, the wide spectrum of anatomical 
connections and the heterogeneous cytoarchitecture are 



related in at least two ways. First, agranular and granular 
regions of the insula establish preferential, though not 
exclusive, connections with cortical regions featuring the 
same architecture, while the dysgranular territory projects 
to regions that range from allocortex to granular isocortex. 
Second, the proximity of different territories in the insula 
dysgranular (Id) with either insula agranular (la), or insula 
granular (Ig) is reflected in their connectivity pattern, in 
such a way that the anterior dysgranular sector has a con- 
nectivity pattern that bears greater resemblance to that of 
la, while the connections of the posterior dysgranular 
insula are more similar to that of Ig [Mesulam and Muf- 
son, 1982b, 1985]. 

Although these studies provided an extensive knowl- 
edge of the anatomical connections of the macaque insula, 
very little is known about the anatomical connectivity of 
the human insula. In this work, we used in vivo probabil- 
istic tractography [Behrens et al., 2003, 2007] on diffusion- 
weighted (DW) magnetic resonance images to investigate 
the connectivity distribution of the insular cortex, where 
we define "connectivity distribution" as the output of 
probabilistic fibre tracking (see Methods, section Probabil- 
istic Tractography for further details). Our aims were to 
analyze the trajectory of connectivity variation across dif- 
ferent insular territories, to compare the results with the 
knowledge derived from tracer-injection studies in maca- 
ques and to evaluate whether the results would suggest a 
relationship between cytoarchitecture and connectivity in 
humans resembling that found in macaques. 

Previous studies have used in vivo probabilistic tractog- 
raphy to successfully recover boundaries between subre- 
gions featuring sharply different connectivity patterns in 
the inferior frontal gyrus [IFG; Anwander et al, 2007; 
Klein et al., 2007], in the lateral premotor cortex [Tomas- 
sini et al., 2007], in the cingulate cortex [Beckmann et al., 
2009], and between supplementary and presupplementary 
motor area (SMA and preSMA, respectively) in the medial 
premotor cortex [mPMC; Johansen-Berg et al., 2004]. In all 
these cases, spectral reordering or hard clustering algo- 
rithms were used, because a clear-cut subdivision was 
expected between subregions with different connectivity 
[Behrens and Johansen-Berg, 2005]. However, in the case 
of the insula, a general consensus about the sharpness of 
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the borders between subregions featuring different connec- 
tivity and cytoarchitecture is still lacking. Brodmann [1909] 
reported two to four main cytoarchitectural subdivisions, 
stressing in several occasions the fact that "there are great 
difficulties in dividing [the insula] into individual fields" 
[Brodmann, 1909: 122]. Other maps of the human insula 
produced in the first half of the 20th century by von Econ- 
omo [1927], Rose [1928], and Brockhaus [1940] reported 
the presence of an anterior agranular and a posterior dys- 
granular sector as well as of a transitional dysgranular sec- 
tor whose location and boundaries vary considerably 
among different maps. More recent investigations in maca- 
ques explicitly reported a gradual transition across differ- 
ent cytoarchitectonic fields [Friedman et al., 1986; Jones 
and Burton, 1976; Mesulam and Mufson, 1982a, 1985] also 
reflected in a gradual differentiation of the intracortical 
myelin [Friedman et al, 1986; Mesulam and Mufson, 
1982a] and associated with a gradual transition of connec- 
tivity [Chikama et al., 1997; Mesulam and Mufson, 1982b, 
1985]. On the other hand, a recent observer-independent 
study on 10 postmortem human brains clearly identified 
boundaries between different cortical types in the insular 
region caudal to the central insular sulcus [Kurth et al., 
2010a]. Finally, at least two studies also reported a similar 
topographical arrangement of different cytoarchitectonic 
fields in humans and macaques [Bonthius et al., 2005; 
Mesulam and Mufson, 1982a]. 

Given the available anatomical evidence, applying a 
hard-clustering technique in the insula would be undesir- 
able, as it would, by definition, split the insula in discrete 
clusters without first assessing whether the data indeed 
suggest the presence of homogeneous clusters separated 
by sudden changes of connectivity. In a previous study 
[Nanetti et al., 2009], we showed that repeated applica- 
tions of fc-means clustering reveal that the location of this 
border is unstable and that a large zone in the middle 
insula cannot be reliably classified using fc-means. There- 
fore, we resort here to the use of a nonlinear dimensional- 
ity reduction technique, Laplacian eigenmaps (LE) [Belkin 
and Niyogi, 2002], to represent insular voxels as points on 
a plane, in such a way that the distance between each pair 
of voxels would quantify the similarity of their respective 
connectivity maps. This technique is similar to the spectral 
reordering algorithm previously used in other connectiv- 
ity-based parcellation studies [Johansen-Berg et al., 2004, 
2008; Klein et al., 2007], but it also takes advantage of the 
information about the nearest neighbors of each connectiv- 
ity feature vector. Considering the spatial arrangement of 
the seed voxels according to the LE of their tractography 
feature vectors allowed us to investigate whether connec- 
tivity varies gradually across the insula or whether it is 
characterized by the presence of sharply delimited 
clusters. 

LE of the insular connectivity maps reveals a main tra- 
jectory of connectivity variation along the rostrocaudal an- 
atomical axis, spanning from the most anterior insula, 
encompassing the anterior short gyrus as well as the ven- 



tral anterior insula around the limen, to the dorsal and 
posterior part of the long gyri. The distribution of the 
brain regions reached by fiber tracking seeded in different 
insular territories is mostly in agreement with tracer-injec- 
tion studies in macaques and suggests a relationship with 
the cytoarchitecture. Finally, LE reveals significant differ- 
ences in the structure of connection variability when com- 
paring insula and mPMC: sharp changes in connectivity 
were detected in the mPMC, while in the case of the 
insula, the connectivity variability was organized in a 
more continuous manner along the identified rostrocaudal 
axis. 



MATERIALS AND METHODS 

Ten healthy adult males (median age, 35; range, 27-50) 
took part in the experiment. Limb dominance was 
assessed with a dutch translation of the Edinburgh Hand- 
edness Inventory [Oldfield, 1971]. All participants were 
right-handed (all scores above 60/100). 

The participants gave informed written consent in ac- 
cordance with the Declaration of Helsinki and with the 
study protocol, which was approved by the Medical Ethi- 
cal Commission (METc) of the University Medical Center 
Groningen (NL). 

Image Acquisition and Preprocessing 

DW data were acquired using a single-shot pulsed gra- 
dient spin echo EPI sequence (TR = 7,500 ms, TE = 67 ms, 
and SENSE factor = 2) on a 3-T MR scanner (Intera, Phi- 
lips Medical Systems, Best, the Netherlands) equipped 
with an eight-channel head coil. For each subject, 15 DW 
images were collected using the "medium" gradient 
scheme provided by the manufacturer, with a maximum 
gradient strength of 22 mT m -1 and a fc-value of 1,000 s 
mm" 2 . In addition, one non-DW volume (b = 0 s mm~ 2 , 
referred in the following as a bO image) was acquired. 
Each volume consisted of 76 transverse slices, acquired 
with a 96 x 96 matrix (FOV = 224 x 224 mm, voxel size = 
2.3 x 2.3 mm, slice thickness 2 mm, no gap; reconstructed 
matrix 128 x 128 leading to a final resolution of 1.75 x 
1.75 x 2 mm 3 ; NEX = 2). The total acquisition time for 
this sequence was 4 min and 14 s. In addition, one high-re- 
solution (1 mm isotropic) 3D anatomical scan (T1W, TFE 
CLEAR, TR = 7,600 ms, TE = 3.6 s, and flip angle = 8°) 
was acquired. 

All the preprocessing stages were performed by means 
of the software provided with FSL 4.1.1 [Smith et al., 2004; 
Woolrich et al., 2009] (available at www.fmrib.ox.ac.uk/ 
fsl). Preprocessing of DW images involved realigning all 
the DW volumes to the bO to correct for movements and 
for eddy currents. Subsequently, the skull was automati- 
cally removed from both the bO and the anatomical scan 
using Brain Extraction Tool [Smith, 2002], and the bO was 
registered to the anatomical scan [affine registration, 7 
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degrees of freedom (DOF)] and to the MNI template 
(affine registration, 12 DOF, target voxel resolution: 2x2 
x 2 mm) by means of FLIRT [Jenkinson et al., 2002]. 

Definition of the Insular Region of Interest 

The insular cortex was drawn on the anatomical image 
realigned to the bO and taken into the corresponding reso- 
lution, using the inverse of the previously calculated ma- 
trix to register the bO to the Tl-weighted image. The 
location of the drawn ROIs was later inspected by overlay- 
ing them onto the bO image. A study from Naidich [Nai- 
dich et al., 2004] showed the reliability of MR images to 
identify anatomical landmarks for the insular region. We 
rely on the landmarks described in that work to draw our 
ROI, using the terminology mediated by Tiire et al. [1999] 
(see Fig. 1). 

The anterior limit of the insula was identified in the an- 
terior peninsular sulcus, at the level of the junction 
between the internal and external capsule, extending ven- 
trally until the subapical insular pole. The posterior limit 
was the most dorsocaudal extent of the anterior long gyrus 
of the insula. The superior periinsular sulcus was defined 
the dorsal boundary of the ROI. Ventrally and caudally, 
the inferior periinsular sulcus (IPS) was drawn in poste- 
rior-to-anterior direction, as far as the junction with the 
temporal operculum (at the level of the limen insulae). No 
anatomical landmark can be used to delimit the insular 
cortex from the adjacent posterior orbitofrontal cortex 
[OFC; Heimer and Van Hoesen, 2006; Mesulam and Muf- 
son, 1982a, 1985]; therefore rostral to the limen insulae, we 
chose the ideal plane connecting the insular pole with the 
most rostral location where the junction with the the tem- 
poral operculum could be identified in order to draw the 
anterior-ventral boundary. Laterally, the insula was delim- 
ited by the CSF, while medially we chose those voxels at 
the interface with the underlying white matter as a bound- 
ary of the seed region. Given the limits imposed by the 
voxel resolution as well as by the morphology of the ven- 
tral insular cortex with respect to the claustrum, it is possi- 
ble that some voxels in the ventral insular ROI may 
include some part of the claustrum (see Discussion: A 
Comment on the Claustrum). 

Probabilistic Tractography 

Probabilistic fiber tracking was performed with the 
FMRIB's diffusion toolbox [Behrens et al, 2003] from 
each insular voxel to the ipsilateral hemisphere using a 
model that accounts for the presence of either one or 
two fibers in each voxel [Behrens et al., 2007]: at each 
step, during the fiber-tracking process, a fiber orientation 
was chosen from each of the estimated fiber populations 
supported by the data, and the pathway was tracked 
along the fiber orientation closest to the one chosen in 
the previous tracking step; pathways looping back on 
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Figure I. 

Definition of the insula ROI. A: Surface rendering of the left 
brain in one subject at approximately mid-way through the cort- 
ical gray matter (cortical layer 4) [Van Essen, 2005]. The coordi- 
nates (0,0) indicate the Y and Z location of the anterior 
commissure (AC). Large parts of the frontal, parietal, and tem- 
poral opercula were removed to show the location of the insula 
and the landmarks chosen to draw the ROI: anterior periinsular 
sulcus (aps), superior periinsular sulcus (sps), and inferior periin- 
sular sulcus (ips). The dotted line connecting the limen insulae 
(li) with the insular pole (ip) represents the ideal boundary 
between the anterior ventral insula and the posterior orbito- 
frontal cortex. A detailed description of the drawing of the 
insula ROI in the subjects' volumes is provided in the main text. 
B: Surface rendering of the insular cortex in the same subject. 
This subject was chosen to show the morphology of the insula, 
because it presents the most common configuration of the 
short and long gyri according to the dissection studies con- 
ducted by Tiire [1 999]. Abbreviations: alg, anterior long gyrus; 
aps, anterior periinsular sulcus; asg, anterior short gyrus; cis, 
central insular sulcus; ia, insular apex; ip, insular pole; ips, infe- 
rior periinsular sulcus; li, limen insulae; msg, middle short gyrus; 
pig, posterior long gyrus; psg, posterior short gyrus; sps, supe- 
rior periinsular sulcus. 
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themselves were terminated. Five thousand samples 
were drawn from the connectivity distribution of each 
seed voxel, yielding a connectivity map recording the 
location of the target voxels reached by probabilistic trac- 
tography as well as the number of samples crossing each 
target voxel. These values were further multiplied by the 
mean length of the pathways connecting each seed with 
each target voxel [Beckmann et al., 2009; Eickhoff et al., 
2010; Tomassini et al., 2007] to ensure that the results of 
tractography, and the following calculation of similarity 
between connectivity maps, were not biased toward 
nearby connections due to the mere anatomical 
proximity. 

It is important to specify how we will frame the results 
of probabilistic tractography in the context of studying the 
anatomical connectivity of the insular cortex. In the follow- 
ing, we will refer to the results of probabilistic tractogra- 
phy as "connectivity maps" and to the samples sent from 
a seed voxel to a target voxel as "projections." The proper 
use of these terms belongs to the realm of animal studies, 
in which tracers can be injected in specific brain regions to 
label axonal pathways, distinguish afferents and efferents, 
and stain single neural termination. In MR tractography, is 
it not possible to achieve this level of resolution; we can 
only infer the location and shape of entire white-matter 
bundles from the effect of the displacement of water mole- 
cules on the MR signal at the voxel level, estimate the 
uncertainty associated with this signal, and finally calcu- 
late a distribution of possible directions to follow during 
tractography. Therefore, when discussing the results of 
probabilistic tractography of the insula, we will use the 
term "projections" to refer to those voxels in the brain 
where samples were recorded after performing probabilis- 
tic tractography from a given seed voxel or seed region, 
the location of these samples in the brain being repre- 
sented by a "connectivity map." In doing so, we comply 
with previous connectivity-based parcellation studies [Beh- 
rens et al., 2003; Croxson et al., 2005; Rushworth et al., 
2006; Tomassini et al., 2007] in which terms such as "con- 
nectivity" and "projections" were used to describe the 
results of probabilistic tractography after strict definitions 
of the terminology. 



LE of the Insula 

Because we aimed at analyzing the structure of the con- 
nectivity variability across insular seed voxels, we first 
needed a quantitative description of the similarity between 
the connectivity map of each and every other seed voxel. 
The similarity between the connectivity maps of two seed 
voxels was quantified in terms of the location of the target 
voxels in the brain reached during probabilistic tractogra- 
phy as well as by the number of samples reaching each 
target voxel. This part of the procedure has been detailed 
in previous works [Anwander et al., 2007; Beckmann 
et al., 2009; Johansen-Berg et al., 2004; Klein et al., 2007; 



Tomassini et al., 2007]. In brief, for every subject and every 
hemisphere, the connectivity map of each seed voxel was 
reshaped in a vector of length equal to the number of vox- 
els in the ipsilateral hemisphere. After reducing the 
dynamic range of the values by logarithmic transformation 
[Anwander et al, 2007], the similarity between each two 
(vectorized) connectivity maps, calculated using the Pear- 
son correlation coefficient, was stored in a square matrix A 
of dimensions n x n, where n is the number of seed voxels 
in the insula of each subject. Each row (or column) of the 
matrix is therefore a feature vector for one seed voxel, rep- 
resenting the similarity between that voxel's connectivity 
map and the connectivity map of every other seed voxel 
(Fig. 2, left). 

The final aim in building this matrix is to group to- 
gether seed voxels with similar connectivity maps. How- 
ever, two connectivity maps can be highly correlated 
because of the connections between their respective seed 
voxel, instead of the connections with brain regions out- 
side the seed mask, which are the most interesting for con- 
nectivity-based parcellation studies. To correct for this 
bias, we added to the matrix A an euclidean distance ma- 
trix of the anatomical distance across the seed voxels 
(scaled in the range 0..1 on a single subject basis). In this 
way, the similarity between two connection probability 
maps is reduced as a function of their anatomical proxim- 
ity. The resulting matrix of the feature vectors was then 
fed into the algorithm of the LE [Belkin and Niyogi, 2002] 
to analyze the trajectory of connection variability across in- 
sular territories. 

LE [Belkin and Niyogi, 2002] was used to recover, if 
present, a structure of connectivity variability across all 
seed voxels' connectivity maps. A technical exposition of 
the implementation of LE for probabilistic tractography 
data is provided in the Supporting Information. In brief, 
the algorithm uses the information about each feature vec- 
tor's nearest neighbours — that is, those vectors with the 
smallest euclidean distance to a given vector — to map the 
feature vectors, which are originally points in R" (a space 
in n dimensions) as points in R 2 (a plane), so that the 
proximity in R 2 reflects the proximity in R". A graph was 
calculated from the matrix A, connecting each two feature 
vectors whose Eluclidean distance was less or equal than 
the distance required for the graph to be connected. The 
Laplacian matrix of this graph was computed, and the 
eigenvectors associated with the first two smallest nonzero 
eigenvalues of the Laplacian matrix (first and second com- 
ponent in Fig. 2, middle) were used to map the connectiv- 
ity feature vectors of each seed voxel as points on a plane. 
If the connectivity across insular regions varies in a contin- 
uous manner between two extremes with highly different 
connectivity patterns, we expected the algorithm to 
arrange the seeds along a line. On the other hand, if the 
insula features subregions with sharply different connec- 
tivity among each other, we expected the algorithm to 
detect large interruptions on the line or segregated clouds 
of points. 
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Figure 2. 



Computation of the Laplacian eigenmap of the insular connectiv- 
ity feature vectors in a single subject. Left column: Similarity ma- 
trix A, containing the correlation values between each and every 
other seed voxel's connectivity maps calculated by probabilistic 
tractography. Colors indicate the value of the Pearson correla- 
tion coefficient R. Middle column: The matrix A was fed into the 
algorithm of the Laplacian eigenmaps (see details in the main 
text and in the Supporting Information) in order to assess the 
presence of a structure of connectivity variability across all 
seeds. The algorithm maps the feature vectors, which are origi- 
nally points in R" (a space in n dimensions) as points in R 2 (a 
plane), so that the proximity in R 2 reflects the proximity in R". 
The recovered one-dimensional structure suggests to the pres- 
ence of a variation between two extremes with very different 
connectivity, while the absence of large discontinuities (further 
quantified in the comparison with the medial premotor region — 
See Methods and Results sections: gradual versus clustered tra- 
jectory of connectivity variability, and Fig. 5) suggests that the 
variation of the connectivity patterns across insular voxels was 



gradual rather than organized in clusters with sharply different 
connectivity patterns. The distance from the minimum value on 
the x axis was calculated for each voxel, normalized in the range 
0.. I and mapped on the anatomical surface of each subject's 
insula (Right column). Note that the distances in the Laplacian 
eigenmap only reflect distances across feature vectors (and ulti- 
mately across the connectivity patterns of each seed voxel), 
while they do not reflect per se distances in millimeters in the 
anatomy: the topographical organization of values assigned to 
the seed voxels emerges only when those distances are mapped 
onto the anatomy. For the purposes of visualization, Caret 
(http://www.nitrc.org/projects/caret/) was used to segment gray 
and white matter and produce fiducial surface representation 
[Van Essen, 2005] of the insula at approximately mid-way 
through the cortical gray matter (cortical layer 4) [Van Essen 
et al., 200 1]. The coordinates (0,0) indicate the Yand Z location 
of the anterior commissure. [Color figure can be viewed in the 
online issue, which is available at wileyonlinelibrary.com.] 



To clarify the subsequent steps of the analysis, we need 
to anticipate that in all subjects, the mapping produced by 
LE revealed the presence of a one-dimensional structure, 
suggesting a trajectory of connection variability between 
two extremes with very different connectivity feature vec- 
tors. Importantly, higher-dimensional LE produced using 
three and four components did not alter the topology of 
the structure. The distance of each point on the line from 
one extreme was therefore calculated and overlaid on the 
anatomy of each subject (Fig. 2, right) to evaluate the topo- 
graphic distribution of the recovered trajectory of connec- 
tion variability. For visualization purposes, these maps 
were then transformed in the space of the standard MNI 
template (affine transformation, 12DOF, and nearest neigh- 
bour interpolation), and the median value for each insular 
voxel belonging to at least five subjects (half of the partici- 



pants) was calculated in order to produce a map of the av- 
erage trajectory of connectivity variability in the insula 
across subjects. By considering only the voxels in MNI 
space that belonged to at least five subjects, we achieved 
the best approximation to the location and boundary of 
the insula as they appear in the ICBM152 brain template 
(http://www.loni.ucla.edu/ICBM/). To evaluate the topo- 
graphic distribution of insular connectivity in relation to 
the results of LE, we considered all the target voxels 
where 50 or more tractography samples (calculated on the 
connectivity maps of each seed voxel for each subject sep- 
arately) had been recorded in at least half of the partici- 
pants. A surface map was then computed, color-matching 
each region of the cortical surface to the insular territory 
from which the largest amount of samples, averaged over 
all the subjects, was recorded during tractography. To do 
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so, the connectivity maps of each subject were grouped in 
100 bins of equal size, according to the trajectory of con- 
nection variability found by LE on a single subject basis. 

Gradual Versus Clustered Trajectory of 
Connectivity Variability 

Tracer injection studies in macaques suggest that the 
insula features a gradual transition of anatomical connec- 
tivity patterns [Chikama et al., 1997; Mesulam and Muf- 
son, 1982b], spanning the rostrocaudal anatomical axis. 
When considering the results of probabilistic tractography, 
LE can provide indications about the relative smoothness 
of the transition among different tractography patterns in 
the insula: for a given seed region, if the transition among 
tractography patterns is relatively smooth, the markers in 
the low-dimensional representation (plane) identifying dif- 
ferent seed voxels (Fig. 2, middle) would lie on a continu- 
ous line with relatively small spatial discontinuities; on the 
other hand, if the examined region of interest is composed 
of different sharply segregated clusters, the structure of 
connectivity variability would display large spatial gaps 
among clusters. To assess which of these two configuration 
would better characterize the tractography results in the 
insula, we compared the LE of the insula with those of the 
mPMC, a region where both animal studies (see Luppino 
et al. [1993] and Picard and Strick [1996] for a review) and 
in vivo connectivity-based panellation studies in humans 
[Anwander et al, 2007; Johansen-Berg et al, 2004; Kim 
et al., 2010] evidenced the presence of sharp changes 
across connectivity patterns. We hypothesized that LE of 
the mPMC would have displayed large gaps on the trajec- 
tory of connectivity variability in correspondence to the 
boundaries between clusters, and these gaps would have 
been quantitatively larger than those present in the LE of 
the insula. 

We calculated LE of the connectivity feature vectors for 
the mPMC region defined in the same way as in previous 
connectivity-based panellation studies [Anwander et al., 
2007; Johansen-Berg et al., 2004]: one single slice on the 
parasagittal plane at X(MNI) = —2, extending from 
Y(MNI) = -22 to y(MNI) = 30 on the rostrocaudal axis 
and from above the cingulate sulcus to the dorsal crest of 
the medial wall of the hemisphere on the dorsoventral 
axis. The analysis was also repeated for the right hemi- 
sphere at X(MNI) = 2. To compare the trajectories of con- 
nection variability of insula and mPMC, we calculated 
the distance of each point in the LE from one extreme, 
and after smoothing the values with a moving average 
(window = five points), we calculated the size of the 
gaps as the distance between the wth point and the 
(n— l)th point along the trajectory of connection variabili- 
ty. To account for the different dimensions of the insula 
and mPMC, the gap values for the SMA region were di- 
vided by the ratio between the size of the insula and the 
size of the SMA/preSMA ROI within each subject. The 



maximum gap was calculated as the median of the high- 
est (1%) gap sizes: this value should inform about the 
presence of transitions between clusters and was there- 
fore expected to be higher for the mPMC. This hypothesis 
was tested separately in each hemisphere with the Wil- 
coxon signed rank test. 

Topographic Distribution of Insular 
Connection Probabilities in Relation to 
the Cytoarchitecture 

In macaques, a strong relationship has been found 
between the cytoarchitecture of different insular regions 
and the cytoarchitecture of the regions in the brain to 
which they are connected [Mesulam and Mufson, 1982b, 
1985; Mufson and Mesulam, 1982]. Using current neuroi- 
maging techniques, it is still difficult to obtain a detailed 
in vivo quantification of the cytoarchitecture of different 
cortical regions (but see Eickhoff et al. [2005a,b], Fatterpe- 
kar et al. [2002], Geyer et al. [2011]). On the other side, 
macroanatomical landmarks do not always match with 
cytoarchitectonic borders [Amunts et al., 1999, 2004, 2005; 
Caspers et al, 2006; Geyer et al., 1999, 2000, 2004; Grefkes 
et al., 2001; Scheperjans et al., 2008]. The standardized 
probabilistic cytoarchitectonic maps available in the SPM 
Anatomy Toolbox [Eickhoff et al., 2005b] currently pro- 
vide the best available estimation of the cytoarchitectonic 
distribution of several cortical regions across subjects. 
Among the currently mapped regions, some are particu- 
larly relevant to the study of insular connectivity, because 
the topographical distribution of their projections in the 
insula is known from animal studies, and because they 
span a wide spectrum of cortical types. We used the 
maps from the cytoarchitectonic atlas of Juelich [Eickhoff 
et al., 2005b] to compare the results of probabilistic trac- 
tography in the human insula with the knowledge about 
insular connectivity derived from tracer injection studies 
in macaques and to assess the relationships between the 
trajectory of connectivity variability recovered in the 
insula with probabilistic tractography, and the trajectory 
of cytoarchitectural variability known from animal 
studies. 

We chose as targets the maximum probability maps 
(MPM, [Eickhoff et al., 2005a,b] of the probabilistic 
cytoarchitectonic maps for BA6 [Geyer, 2004], SI (Area 1, 
2, 3a, 3b) [Geyer et al, 2000; Grefkes et al., 2001], BA44 
and BA45 [Amunts et al., 1999], inferior parietal lobule 
(IPL: PF, PFm, PFcm, PFt, PGp, Pfop, and PGa; [Caspers 
et al., 2006]), superior parietal lobule (SPL: 5M, 5Ci, 5L, 
7M, 7PC, 7A, and 7P; [Scheperjans et al, 2008]), and a vol- 
ume, hereafter referred to as Limbic ROI comprising hip- 
pocampus (CA, FD, HATA, and SUB), amygdala (CM, SF, 
and LB), and enthorinal cortex (EC) [Amunts et al., 2005]. 
The MPM of these regions was transformed in the single- 
subject DW space using the inverse of the previously cal- 
culated transformation matrix taking the single-subject bO 
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image in the MNI space (affine transformation, 12DOF, 
and nearest neighbour interpolation), and the connection 
probability of each seed voxel to each target region was 
calculated as the proportion of the samples from the con- 
nectivity distribution of that seed, which was sent to all 
targets. For each target region, voxels in the insula with a 
connection probability lower than 0.05 were excluded 
from further analyses. For the purposes of visualization 
and qualitative estimation, the resulting seven connection 
probability maps for each subject were transformed into 
the standard MNI space (affine transformation, 12DOF, 
and trilinear interpolation), and the average (median) con- 
nection probability maps across subjects were calculated in 
all insular voxels belonging to at least half of the partici- 
pants (five subjects — for the reasons previously detailed in 
Methods section: LE of the insula). 

To investigate potential asymmetries in insular con- 
nections between right and left hemisphere, we consid- 
ered separately the amount of samples sent from the 
insular cortex to each target, and the extent of insular 
cortex found to be connected to each target region. For 
each subject, voxels sending less than 50 samples to 
each target were discarded to remove potential biases of 
the results due to seed voxels with a low probability to 
be connected to the target. Differences between hemi- 
spheres in the average number of samples and in the 
extent of the connected insular region to each target 
were computed by means of paired t tests and corrected 
for multiple comparison with a false discovery rate of 
0.05. To investigate the spatial variability of insular con- 
nections across subjects and between hemispheres, we 
computed maps showing for each insular voxel the 
number of subjects who displayed an above-threshold 
number of tractography samples to each of the exam- 
ined targets and plotted the associated centers of gravity 
(COG). For each target, we also used Levene and 
Brown-Forsythe tests to compare the variability in the 
distance of every subject's COG from the mean COG 
between the two hemispheres. 

RESULTS 
LE of the Insular Cortex 

In each subject and hemisphere, Laplacian eigenmaps 
(LE) identified a topographically organized trajectory of 
connectivity along the rostrocaudal axis. One extreme of 
this trajectory was located in the anterior insula, in a terri- 
tory including the ventral anterior insula around the limen 
as well as the anterior short gyrus; the other extreme was 
located in the dorsal and caudal extent of the long 
gyri and in the adjacent part of the central insular sulcus 
(Figs. 3 and 4). In the representation of LE, these extremes 
are therefore the two insular fields differing most in terms 
of their respective connectivity patterns. 

To evaluate the topographic distribution of the connec- 
tivity patterns from different insular territories, we 



masked all the regions on the cortical surface reached by 
probabilistic tractography in at least half of the partici- 
pants (five subjects), and we color-matched each location 
on the surface with the insular location from which the 
largest amount of samples, averaged over all the sub- 
jects, were recorded during tractography (Fig. 4, bottom 
row). The cortical regions receiving most of the projec- 
tions from the anterior insula included the OFC, pars 
orbitalis and several locations in pars triangularis of the 
IFG, and the dorsal part of the temporal pole. Further 
analyses (see below: Topographical and cytoarchitectonic 
specificity of connectivity patterns across insular territo- 
ries) found the centromedial nucleus and laterobasal 
complex of the amygdala to be connected prevalently 
with the anterior insula and specifically with the most 
ventral and anterior insular territory around the limen. 
On the other hand, most of the tractography samples 
from the posterior dorsal insula reached the parietal 
lobe, including SI, SII, and the posterior parietal lobule 
as well as posterior regions of the temporal lobe in the 
superior temporal gyrus, the superior temporal sulcus, 
and the middle temporal gyrus. Samples from the 
insula, and especially from this posterior portion, were 
also recorded in extrastriate regions of the occipital cor- 
tex both medially (lingual gyrus) and laterally (putative 
BA19; about this see Discussion section: A comment on 
the claustrum). Regions in the frontal lobe that were 
mostly connected with the posterior dorsal insula 
included cortical territories in the premotor cortex, in 
pars triangularis and pars opercularis of the IFG as well 
as in the left dorsolateral prefrontal cortex (putative 
BA46). Middle insular territories displayed a mixed con- 
nectivity pattern between anterior and posterior insula. 
Cortical regions receiving the largest amount of samples 
from middle insular territories were located mostly in 
the lateral and medial premotorcortex as well as in 
regions of the temporal and parietal lobe adjacent to 
those reached by the posterior dorsal insula. Several ter- 
ritories in the IFG and OFC were also mostly reached by 
samples from the middle insula. 

Gradual Versus Clustered Trajectory of 
Connectivity Variability 

We analyzed the LE of the insula and of the mPMC to 
test the hypothesis of a smoother transition across connec- 
tivity patterns in the insula compared to a region like 
mPMC where previous studies found clusters of seed vox- 
els featuring sharply different connectivity patterns. Wil- 
coxon signed rank test revealed that the size of the gaps 
along the trajectory of connectivity variability defined by 
the LE was significantly higher for the mPMC than for the 
insula in both hemispheres, suggesting the presence of 
comparatively sharper clusters in the examined medial 
frontal ROI with respect to the insula: right mPMC median 
gap = 0.0198, right insula median gap = 0.0125: W mPMC = 
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Laplacian eigenmaps of the insula for each subject: Laplacian 
eigenmaps of the connectivity feature vectors derived from the 
in vivo probabilistic tractography in each subject and hemi- 
sphere. For purposes of visualization, Caret (http://www.nitr- 
c.org/projects/caret/) was used to segment gray and white 
matter and produce fiducial surface representation [Van Essen, 
2005] of the insula at approximately mid-way through the corti- 
cal gray matter (cortical layer 4) [Van Essen et al., 200 1]. The 
final images of the surfaces and the mapping of the metrics 



derived from the Laplacian eigenmaps were rendered by means 
of Paraview (http://www.paraview.org/). Values on the surface 
were mapped from the volume containing the results of Lapla- 
cian eigenmaps and quantifying the distance along the main tra- 
jectory of connectivity variability: each vertex on the surface 
was assigned the value of the closest voxel in the volume. 
[Color figure can be viewed in the online issue, which is avail- 
able at wileyonlinelibrary.com.] 
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138, Z = 2.46, P < 0.014.; left mPMC median gap = 0.0218, 
left insula median gap = 0.0057: W mPMC = 149, Z = 3.29, 
and P < 0.0011 (see Fig. 5). 

Topographical and Cytoarchitectonic Specificity 
of Connectivity Patterns Across Insular 
Territories and Between Hemispheres 

Tracer-injection studies in macaques evidenced a topo- 
graphic similarity between the trajectory of variation in 
the cytoarchitecture, organized radially around the ventral 
anterior agranular insula, and the trajectory of variation of 
insular projections, organized along a rostrocaudal axis 
[Chikama et al, 1997; Mesulam and Mufson, 1982b]). The 
rostrocaudal gradient of connectivity variability recovered 
in our study with LE was therefore expected to reflect 
changes in the connection probability with regions featur- 
ing different cortical types. To assess this hypothesis, we 
used 7 MPM of cytoarchitectonically defined regions in 



the atlas of Juelich [Eickhoff et al., 2005b] spanning a wide 
spectrum of cortical types ranging from agranular (amyg- 
dala, hippocampus, EC, and BA6) to dysgranular (BA44) 
to granular cortex (BA45, SI, SPL, and IPL), and we quan- 
tified the connection probability, averaged across subjects, 
of each voxel in the insula with each MPM. For example, 
the connection probability of a voxel in the insula with 
BA44 is quantified as the number of samples from that 
voxel reaching BA44, divided by the total number of sam- 
ples reaching all the seven explored target regions. This 
measures the selectivity of the connections from insular 
voxels to targets, rather than the absolute number of sam- 
ples (the latter was quantified previously, and shown in 
Fig. 4). Results of this analysis are shown in Figure 6 A, 
overlaid onto the insular surfaces of the single-subject 
MM template. 

The connection probability maps show that some insular 
territories are characterized by relatively exclusive connec- 
tions with specific targets, while others send a comparable 
number of seeds to all examined targets. In particular, the 



Figure 4. 



Average trajectory of connectivity variability and maximum num- 
ber of tractography samples per target voxels. Average trajec- 
tory of connectivity variability (upper row): For visualization 
purposes, after calculating the Laplacian eigenmaps of the insular 
connectivity for each subject and each hemisphere, and replot- 
ting them into each subject's anatomical space, the volumes 
were transformed into the standard MNI space (see details in 
Methods section: Laplacian eigenmaps of the insula), and the 
median value for each insular voxel belonging to at least half of 
the participants was calculated and mapped onto the insular sur- 
face of the MNI single subject template [Holmes et al., 1998]. 
This illustration serves to show the rostrocaudal topographical 
arrangement of the recovered trajectory of connectivity variabil- 
ity, spanning from the most anterior insular territory, in the an- 
terior short gyrus and in the ventral anterior insula around the 
limen, to the posterior insular territory, in the dorsal caudal 
long gyri and adjacent central insular sulcus. The same trajec- 
tory of variability was recovered in all subjects, as it can be 



appreciated from Figure 3. Additional analyses, comparing the 
single subject's Laplacian eigenmaps of insula and medial premo- 
tor cortex, were performed to test for the gradual versus clus- 
tered organization of connectivity patterns (see Methods and 
Results sections: Gradual versus clustered trajectory of connec- 
tion variability; Discussion section: Gradual variation of connec- 
tivity patterns). Maximum number of tractography samples per 
target voxels (lower row): We thresholded the single subject's 
connectivity map for each seed voxel to a minimum of 50 sam- 
ples, and we considered only target brain voxels reached in at 
least half of the participants (five subjects). Then, we color- 
matched each of the surviving target brain voxels with the insu- 
lar location, derived from the average trajectory of connectivity 
map, from which it received the maximum number of samples, 
averaged over all subjects present in that location. [Color figure 
can be viewed in the online issue, which is available at 
wileyonlinelibrary.com.] 



Figure 5. 



Comparison between insula and medial premotor cortex. A: 
Laplacian eigenmaps identified the main trajectory of connectiv- 
ity variability in the insula along the rostrocaudal axis. The ab- 
sence of large gaps in the Laplacian eigenmap suggests the 
presence of a smooth transition across different connectivity 
patterns. B: On the other hand, the medial premotor cortex 
(mPMC) features a sharp transition between clusters with differ- 
ent connectivity patterns, which are evidenced by the presence 
of large gaps in the corresponding Laplacian eigenmap. The Lap- 
lacian eigenmap of the mPMC is overlaid onto the MNI single 
subject template to show that the transition is located on the 
plane through the anterior commissure (indicated by the cross- 



hairs). The size of the largest gaps in the Laplacian eigenmaps of 
insula and mPMC was therefore taken as a measure of the pres- 
ence of sharp transition in connectivity patterns and compared 
between the two seed regions across all subjects. C: Boxplots 
(median, interquartile range, and most extreme values within 1.5 
times the interquartile range) of the maximum detected gap in 
the Laplacian eigenmaps of insula and mPMC in both hemi- 
spheres. Wilcoxon signed-rank test showed a significant differ- 
ence between the medians of the maximum detected gap in the 
insula and the mPMC. [Color figure can be viewed in the online 
issue, which is available at wileyonlinelibrary.com.] 



♦ 2015 ♦ 



♦ Cerliani et al. ♦ 



Limbic y" 



Limbic 



BA44 



o ■ 



BA45 
1 



^^^^ 



BA44 




BA45 



BA6 



0.2 I 



SI 

0.4 I 



BA6 



0.1 1 




SI 



0.61 



SPL 



0.4 



IPL 



seed/target connection probability 



SPL 



0.8 



IPL 



B1 




>."Ji ■ i " . i ■ 1 1 ;■■ 



B2 



B3 



I - 




Figure 6. 



♦ 2016 ♦ 



♦ Rostrocaudal Connectivity Variations in the Insula ♦ 



ventralmost anterior insula around the limen, location of 
the putative agranular insula, sends the vast majority of its 
samples to the limbic ROI, the dorsal anterior insula to 
BA45, and the caudal portion of the long insular gyri to 
the posterior parietal lobule. The middle insula, on the 
other hand, sends comparable amounts of samples to all 
but the limbic target. 

When comparing the number of samples sent from the 
insular cortex to each of the examined targets in the RH 
and LH hemisphere, we found significant differences in 
BA44 (RH > LH, t = 4.33, P < 0.002) and SPL (RH > LH, f 
= 3.24, P < 0.011). On the other hand, the extent of insular 
cortex found to be connected to each target significantly dif- 
fered for the Limbic target (LH > RH, t = 3.41, P < 0.0077) 
and for BA6 (LH > RH, t = 2.543, P < 0.032; Supporting In- 
formation Fig. S3). The latter result was also reflected in the 
maps showing for each insular voxel the number of subjects 
with an above-threshold number of samples to that target 
(Fig. 7, up). It is noteworthy the fact that in these two 
regions as well as in SI, there is a considerable spread of 
the centers of gravity of the connected insular region in the 
right hemisphere, suggesting a degree of variability across 
subjects in the extent of the connected insular cortex, even 
though the amount of samples sent to the respective target 
turned out to be comparable. However, we note that no sig- 
nificant difference was found when comparing between the 
two hemispheres, and for each target region, the variability 
in the distance of every subject's COG from the mean COG 
(Levene test: all P > 0.184; Brown-Forsythe test: all P > 
0.174). Finally, in the case of SI, IPL, and left BA45, Figure 7 
(up) suggests that across subjects, a consistent amount of 
tractography samples were sent from two spatially different 
insular territories. 



DISCUSSION 

We examined the anatomical distribution of insular con- 
nections in humans by means of in vivo probabilistic trac- 
tography in order to compare the results with the 
evidence from tracer-injection studies in macaques. LE on 
the results of probabilistic tractography showed that the 
main trajectory of connectivity variation in the human 
insula is organized along the rostrocaudal axis. The subse- 
quent connectivity-based classification analysis suggested 
a relationship between the topographical variation of con- 
nectivity patterns and the variation of cytoarchitecture, 
similar to what has been found in macaques. The compari- 
son between the LE of insula and mPMC suggests that the 
regional transition among different tractography patterns 
in the insula is smoother than in other brain regions. 

In the following, we will first discuss these main find- 
ings. Also, because patterns of anatomical connectivity are 
important indicators of the functionality of different brain 
regions, we will discuss our results in comparison with 
recent findings about the involvement of the insula in 
speech production and in the experience of pain. In the 
speech domain, we will assess the specificity of projections 
from the anterior insular cortex to brain regions involved 
in speech production. In the pain domain, we will compare 
the locations in the insula where a prevalence of either lim- 
bic or somatosensory connections was found by probabilis- 
tic tractography, with the location of insular territories 
specialized for different aspects of the pain experience. 
Finally, we will discuss our results in relation with the 
findings from recent meta-analyses on the functional orga- 
nization of the insula [Kurth et al., 2010b; Mutschler et al., 
2009], and we will review two recent parcellation studies 
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Connection probability maps with cytoarchitectonically defined 
regions and results from previous fMRI meta-analyses. A: Average 
(median) connection probability for each insular voxel present in 
at least half of the participants (five subjects), calculated for the 
maximum probability maps (MPM) of several cytoarchitectonically 
defined regions in the atlas of Juelich [Eickhoff et al., 2005b]. The 
values are plotted on the insular surface of the MNI single subject 
template. The Limbic ROI (amygdala, EC, and hippocampus), 
BA44, BA45, BA6, primary somatosensory cortex (SI), superior 
(SPL), and inferior (IPL) parietal lobule MPMs were chosen to eval- 
uate the different connection probability of different insular terri- 
tories because of their known connectivity with the insula from 
animal studies as well as for the wide range of cortical types that 
they span. The anterior ventral insula, putative location of the 
agranular insula, had the highest connection probability with the 
limbic ROI, while the highest connection probability with poste- 
rior parietal regions was found in the posterior gyri of the insula, 
putative location of the granular, and adjacent dysgranular insula. 
These two regions featured the two extremes of the trajectory of 



connection probability found with Laplacian eigenmaps (see Figs. 3 
and 4). BI-B3: Maps showing the results of the meta-analysis of 
the functional organization of the anterior insula performed by 
Mutschler et al. ([2009] — reproduced with permission of Elsevier) 
on 58 fMRI studies, reporting in total 159 peaks of activation. Bl 
refers to studies related to peripheral and autonomic changes as 
well as co-activation of voxels in the insula and in the amygdala; B2 
refers to studies on auditory and language processing; B3 refers to 
studies involving hand movements. B4: Result of the meta-analysis 
performed by [Kurth et al., 2010] on coordinates from 1768 stud- 
ies retrieved from the BrainMap database [Fox et al., 2005a; Laird 
et al., 2005, 2009] and Pubmed. In the result shown here (adapted 
from Fig. 7 in Kurth et al. [2010] and reprinted with permission of 
Springer), studies from 1 3 functional categories were grouped in 
four domains to investigate the presence of domain-specific func- 
tional specialization in the entire insular cortex: sensorimotor 
(red), cognitive (green), chemical sensory (yellow), and social 
emotional (blue). [Color figure can be viewed in the online issue, 
which is available atwileyonlinelibrary.com.] 
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on the insula based on resting-state functional connectivity 
[Cauda et al., 2011; Deen et al., 2010]. 



Trajectory of Connectivity Variability in the 
Insular Cortex 

Our results are in line with the evidence from tracer- 
injection studies in macaques in several respects. First, LE 
revealed a relatively smooth transition across connectivity 
patterns, topographically arranged along the rostrocaudal 
anatomical axis. Several cytoarchitectonic studies in maca- 
ques explicitly report that a gradual transition across corti- 
cal types [Friedman et al., 1986; Jones and Burton, 1976; 
Mesulam and Mufson, 1982a, 1985], as well as in the 
amount of intracortical myelin [Friedman et al., 1986; 
Mesulam and Mufson, 1982a], is reflected in the transition 
among different connectivity patterns [Chikama et al., 
1997; Mesulam and Mufson, 1982b, 1985]. Also, tracer- 
injection studies in macaques suggested the presence of a 
gradual change in the insular connections with the amyg- 
dala [Amaral and Price, 1984], the striatum [Chikama 
et al., 1997; Fudge et al., 2005] and the secondary somato- 
sensory cortex [Friedman et al., 1986]. 

Second, the comparison between the topographical 
arrangement of seed and target voxels suggests the pres- 
ence in humans of a relationship between the variation of 
the cortical types and the variation in the anatomical con- 
nections as previously found in macaques. Tracer-injection 
studies evidenced that the la and Ig have very different 
connectivity, often with brain regions featuring the same 
cytoarchitecture [Amaral and Price, 1984; Augustine, 1985, 
1996 for a review; Fudge et al., 2005; Mesulam and Mufson, 
1985; Mufson and Mesulam, 1982]. In our study, as we 
matched each location on the cortical surface, reached by 
probabilistic tractography in at least half of the partici- 
pants, with the insular territories from which they received 
the maximum projections, we observed that the OFC, the 
anterior parts of the IFG as well as the amygdaloid com- 
plex were mostly reached by samples from the connectivity 
distribution of the anterior insular territory encompassing 
the ventral anterior insula around the limen and the ante- 



rior short gyrus, where the la and the adjacent Id were, 
respectively, located in previous histological studies [Bon- 
thius et al., 2005; Mesulam and Mufson, 1982a; Ongur 
et al., 2003]. On the other hand, parietal and posterior tem- 
poral regions, several locations in the lateral premotor cor- 
tex, and in posterior territories of the IFG were mostly 
reached by samples from the posterior insular territory 
encompassing the dorsal part of the long insular gyri and 
the adjacent central insular sulcus, in the location of Ig and 
adjacent Id [Bonthius et al., 2005; Kurth et al., 2010a; Mesu- 
lam and Mufson, 1982a]. 

Third, in macaques, the dysgranular insula projects to 
brain regions ranging from allocortex to homotypical iso- 
cortex [Amaral and Price, 1984; Augustine, 1985, 1996 for 
a review; Fudge et al., 2005; Mesulam and Mufson, 1985; 
Mufson and Mesulam, 1982], and different territories of Id 
display different connectivity patterns according to the 
proximity with either la or Ig [Mesulam and Mufson, 
1985]. In our study, we found that territories lying in the 
middle insula displayed a wide spectrum of connections 
including cortical regions in the parietal and temporal 
lobe, the whole IFG, and scattered locations in the OFC as 
well as widespread connections to the premotor cortex. 
The anterior extreme of the connectivity variability encom- 
passed the putative locations of both la and adjacent Id, 
which were found to have very similar connectivity pat- 
terns in macaques (see Fig. 9 in Mesulam and Mufson 
[1982b]), while samples from more caudal insular territo- 
ries, in the putative posterior Id, were mostly found in 
regions of the cortical surface adjacent to those reached by 
the other extreme of connectivity variability in the most 
dorsal and caudal insula. 



Comparison with Evidence from Studies 
in Macaques 

A relationship between connectivity and cytoarchitec- 
tonic variation was also suggested by the analysis of the 
insular connection probability to 7 MPM taken from the 
probabilistic cytoarchitectonic atlas of Juelich [Eickhoff 
et al., 2005b] and chosen among the available maps in 



Figure 7. 



Lateralization and consistency maps. Up: Consistency maps show- 
ing for each insular voxel the amount of subjects which showed an 
above-threshold amount of tractography samples to each of the 
seven examined targets: Limbic ROI (amygdala, enthorinal cortex, 
and hippocampus), BA44, BA45, and BA6, primary somatosensory 
cortex (SI), superior (SPL), and inferior (IPL) parietal lobule. The 
parasagittal slice showed here was chosen to display the coordi- 
nate corresponding to the maximum overlap across subjects. For 
clarity, only voxels with three or more overlapping subjects are 
displayed. In the case of SI, IPL, and BA45, consistent projections 
across subjects were found in two spatially different insular terri- 



tories (in one or both hemispheres). Black diamonds indicate the 
locations of the center of gravity (COG) of the connected cluster 
in different subjects. Down: Sulcal variability map. After registra- 
tion of each-subject skull-stripped anatomical (T I -weighted) image 
to the MNI single subject template [Holmes et al., 1998], sulcal 
outlines were drawn manually on two single sagittal slices. The 
consistency map of the location of the sulci across is presented 
here overlaid onto the average anatomy across all the participants. 
The outline of each subject's sulci can be seen in Supporting Infor- 
mation Figure S4. [Color figure can be viewed in the online issue, 
which is available atwileyonlinelibrary.com.] 
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order to include regions known to be connected with the 
insula and to span a wide spectrum of cortical types: 
BA44, BA45, BA6, SI, inferior and superior parietal lobule, 
and a limbic ROI, encompassing amygdala, EC, and hip- 
pocampus. The limbic ROI had the strongest connection 
probability to the anterior ventral insula, while SI and IPL 
had the strongest probability of connection with the poste- 
rior dorsal insula. Territories in the putative Id shared 
similar connection probability with all the examined maps 
and were also the target of the maximum connection prob- 
ability for BA44, BA6, and BA45. 

The topographic arrangement of connection probabilities 
in the insular cortex with these regions was also mostly 
consistent with studies in macaques and will be consid- 
ered separately in what follows. 



Limbic connectivity 

Anatomical studies have consistently found that the in- 
sular cortex has extensive limbic and paralimbic connec- 
tions concentrated mostly in its anteroventral territories 
(see Augustine [1985] and Mesulam and Mufson [1985] for 
a review), which have themselves been included in the 
definition of the limbic lobe [Heimer and Van Hoesen, 
2006; Morgane and Mokler, 2006]. To quantify insular con- 
nectivity with other structures of the limbic system, we 
used a cytoarchitectonically defined ROI encompassing 
amygdala, EC, and hippocampus [Amunts et al, 2005]. 
We found the highest probability of connection with these 
regions in the most rostral ventral insular field, corre- 
sponding to the putative agranular insula, in line with ani- 
mal findings [Aggleton et al., 1980; Amaral and Price, 
1984; Kerr et al., 2007; Mohedano-Moriano et al, 2007; 
Mufson et al., 1981; Stefanacci and Amaral, 2002]. 



IFG and premotor cortex 

Comparing the projections to the inferior frontal gyrus 
(IFG) and adjacent ventral premotor cortex in our study 
with evidence from animal studies is complicated by the 
lack of a definitive agreement regarding the homologies 
between the subdivisions of this region in macaques and 
humans [Gerbella et al., 2007; Geyer et al., 2000; Petrides 
and Pandya, 2002; Petrides, 2005; Rizzolatti et al., 1998]. 
Although the exact border separating the macaque homo- 
logues of human BA45, 44, and 6 remains disputed, a cer- 
tain consensus is emerging: the anterior bank of the 
inferior ramus of the arcuate sulcus and the adjacent ros- 
tral convexity are thought to be homologous to human 
BA45, while its posterior bank and adjacent caudal con- 
vexity fall within the homologues to human BA44 and lat- 
eral BA6 [Belmalih et al., 2007; Brodmann, 1909; Geyer 
et al., 2000; Matelli et al, 1985; Nelissen et al., 2005; Pet- 
rides, 2005; Rizzolatti and Arbib, 1998; Vogt and Vogt, 
1919]. 



BA44 and ventral BA6. Mesulam and Mufson [1982b] and 
Mufson and Mesulam [1982] found the caudal bank of the 
ventral arcuate sulcus to be connected with the middle 
insula. Matelli et al. [Matelli et al., 1986] replicated these 
results, also finding the caudal convexity to be connected 
with the middle insula, but in the latter study, the staining 
extended slightly more rostrally in the dysgranular insular 
territory. Consistent with these evidences, we found that 
BA44 and ventral BA6 had their maximum connection 
probability to the dorsal part of the middle insula in the 
territory of the middle and posterior short gyri (see Fig. 6). 
We also note that most of the projections to ventral BA6 
originated from the middle insula (see Fig. 4). 

BA45. Previous studies showed evidence for moderate an- 
atomical connections between dysgranular and granular 
insula and the macaque homologue of BA45 [Gerbella 
et al., 2010; Mesulam and Mufson, 1982b; Petrides, 2006; 
Petrides and Pandya, 2002]. In agreement with these find- 
ings, Figure 4 shows that sectors of the left BA45 receive 
most of the tractography samples from territories in the 
middle and posterior insula. We also observed that the an- 
terior dorsal insula (Fig. 6A), found to be involved in lan- 
guage processing (see Fig. 6B2 and Discussion: 
Comparison with functional studies), sends more samples 
to BA45 than to the other six ROIs we considered for 
quantifying connection probability. This latter finding 
should however be considered in the light of the fact that 
more posterior territories of the insula total substantially 
more samples to BA45 than this anterior region (Fig. 4 and 
Supporting Information Fig. SI). 

Dorsal and medial BA6. In macaques, the posterior insula 
was shown to be connected with both F3 and F6 on the 
medial wall, with the most consistent connections between 
the granular insula and F3 [Luppino et al., 1993]. In our 
study, we found tractography samples from the insula in 
the most dorsal extent of the premotor cortex as well as in 
the supplementary motor area (the human homologous of 
macaque area F3 [Picard and Strick, 1996]). The seed 
regions were mostly located in the putative posterior dys- 
granular and granular insula (see Fig. 4). 

Primary somatosensory cortex 

Tracer injection studies in macaques have shown that SI 
is strongly reciprocally connected with the insula and 
especially with posterior territories spanning the granular 
and adjacent dysgranular fields [Darian-Smith et al., 1993; 
Friedman et al., 1986; Mesulam and Mufson, 1982b]. Also, 
a large portion of Ig was shown to be modality-specific for 
tactile stimuli [Schneider et al., 1993]. Mesulam and Muf- 
son [1985] indicated a bifocal projections to the face and 
leg/foot representation territories of SI in the dorsal Id 
and Ig, possibly underlying a somatotopic organization. 
We found a similar connectivity pattern with projections 
both to the foot and the dorsal crest of the anterior 
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postcentral gyrus, mostly originating in the middle and 
caudal insular cortex. 

Inferior and superior parietal lobule 

Several anatomical works in the last decade have consis- 
tently reported the presence of anatomical connections 
between the posterior parietal cortex and the insula, with 
projections to the IPL, and especially from its anterior por- 
tion, concentrated mainly in the granular insula and to a 
lesser extent in the adjacent dysgranular fields [Andersen 
et al., 1990; Cavada and Goldman-Rakic, 1989; Mesulam 
and Mufson, 1982b, 1985; Neal et al, 1987; Rozzi et al, 
2006]. In our study, we found the posterior dorsal part of 
the insula, in the region of the long insular gyri, to have 
the maximum connection probability with IPL. The maxi- 
mum connection probability with SPL, on the other side, 
was recorded in a region corresponding to the posterior 
long gyrus, in line with the observations from Mesulam 
and Mufson [1982b]. 

Differences between human and macaque insula 

There are currently two strong sources of information 
about the connectivity of the human brain: (1) in vivo trac- 
tography studies of the human brain, with their technical 
shortcoming — that we discuss in the Conclusion section of 
this work — and (2) invasive tracer injection studies in the 
macaque brain, with their problems associated with differ- 
ences between species. Although these two methods are 
each affected by their own shortcoming, the source of 
errors in the two methods is different and independent. 
We and others before us [Anwander et al., 2007; Beck- 
mann et al., 2009; Croxson et al, 2005; Johansen-Berg 
et al, 2004, 2008; Klein et al., 2007; Tomassini et al, 2007] 
therefore seek to reduce the overall uncertainty by generat- 
ing our hypotheses for tractography studies and validate 
our results using the evidence gathered from tracer injec- 
tion studies in macaques. We therefore feel that combining 
two methods with independent problems, namely human 
tractography and macaque tracer studies, remains to date 
the best approach to shed light on human connectivity. 
Until more precise methods will be developed, which can 
deliver the same accuracy of tracer-injection studies, we 
can only highlight the current evidence suggesting actual 
differences across species, and reason about the possible 
implications for tractography results. 

One difference between human and macaque insula 
regards the Von Economo neurons (VEN) [Nimchinsky 
et al., 1999; von Economo and Koskinas, 1925], which are 
present in the insula and in the anterior cingulate cortex 
(ACC) of humans and great apes, while they are absent in 
macaque. These neurons were found in humans especially 
in the agranular insula and gradually decreasing toward 
the superior anterior insula [Allman et al., 2010]. It has 
been hypothesized that these neurons could represent the 



substrate for fast relaying interoceptive, emotional, cogni- 
tive, and sensory information to other parts of the brain, 
particularly for the purpose of generating fast adequate 
response in fast-paced behaviours, such as during social 
interaction [Allman et al., 2010; Craig, 2009]. Indeed, these 
neurons are greatly reduced in pathological conditions 
presenting marked deficits in social cognition such as fron- 
totemporal dementia [Seeley, 2010; Seeley et al., 2006] and 
agenesis of the corpus callosum [Kaufman et al., 2008]. 
Craig proposed that VENs could be involved in the fast 
interconnections between the anterior insula and the ACC 
[Craig, 2009]. A first result obtained with diffusion tensor 
tractography in a 27-year-old male gorilla indicated con- 
nectivity of the insular region cointaining VENs with the 
frontal pole, as well as other parts of frontal and insular 
cortex, the septum and the amygdala, and close examina- 
tion of the reported images also suggests a connection 
with ACC [Allman et al., 2010]. Using very high-spatial re- 
solution DW data (90 micra cubic voxels) on formalde- 
hyde-fixed specimens, researchers from the same lab were 
able to reconstruct connectivity maps between the insula 
and ACC in mouse lemur, macaques, and gorillas (Park 
et al., personal communication). In a recent in vivo tractog- 
raphy study in humans, connections between the anterior 
insula and ACC have been reported [van den Heuvel 
et al., 2009]; however, in that study, the insular and ACC 
ROIs, identified by resting-state functional connectivity 
analysis, apparently extended to adjacent territories (in the 
opercula and in the medial motor wall, respectively — cfr. 
Fig. 6 in van den Heuvel et al. [2009]) and were also 
dilated to a maximum of 4 mm before performing tractog- 
raphy. Another study using probabilistic tractography and 
focusing on the connection of the entire cingulate cortex 
did not report any connections with the insula [Beckmann 
et al., 2009], and the same result was found in the present 
work. We believe that these diverging results are due to a 
limitation in the current methods in tracing fibers through 
major fiber bundles like the corona radiata and the corpus 
callosum, and more generally in regions of high structural 
complexity [Beckmann et al., 2009], as well as from the 
limitations of the quality of the DW data, which are cur- 
rently obtainable with in vivo studies in humans. With 
respect to our study, because our hypothesis was derived 
from evidences in macaque, in which VENs appear not be 
present, it is legitimate to ask in which way this difference 
would affect the trajectory of connectivity variability in the 
insula. Preliminary tractography results from our lab sug- 
gest that LE recovers an analogous difference in the rate of 
change of connectivity patterns in humans and macaques 
when comparing insula and mPMC (see Supporting Infor- 
mation Fig. S2 and [Cerliani, 2009]). Nevertherless, current 
MR tractography techniques can only analyze the connec- 
tivity of large ensembles of neurons and axons and there- 
fore cannot specifically examine the connectivity of VENs. 
Also, to our knowledge, the connectivity of VENs still 
needs to be analyzed by means of traditional tracer injec- 
tion techniques [Allman et al., 2010]. Finally, the 
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connectivity between insula and ACC in humans needs to 
be clarified by further analyses, given the divergence of 
the results in the current literature. It is therefore prema- 
ture at this stage to attempt to answer to the question of 
whether and how the presence of VENs in humans but 
not in macaque would have affected the trajectory of con- 
nectivity variability, which can be recovered with probabil- 
istic tractography. 

Another notable difference between the functional neuro- 
anatomy of human and macaque insula was found by 
studying taste representation in the insula. In macaques, 
anterograde tracer injections have shown that the primary 
gustatory cortex is located in the most rostral part of the 
superior limiting sulcus and more precisely in a small 
granular region at the insular /opercular junction of the 
frontal lobe [Pritchard and Norgren, 2004; Pritchard et al., 
1986; Small, 2010]. On the other side, several functional 
neuroimaging studies in humans (reviewed in Small et al. 
[1999]) showed that taste response is localized in other, 
more posterior regions of the insula /frontal operculum, 
and that the most anterior dorsal part of the insula may be 
involved in a more integrative activity in the context of 
feeding behaviour [Small, 2010; Veldhuizen et al, 2008]: in 
the context of taste studies, Jabbi and colleagues [Jabbi 
et al., 2007] showed that the anterior dorsal insula is 
involved in both experience of disgusting tastants and ob- 
servation of disgusted facial expressions in others and that 
the different sensitivity of the participants to the feelings of 
others predicted the activity of the same insular region 
when observing both pleased and disgusted facial expres- 
sions of others. Rather than a caudal "shift" of the primary 
gustatory cortex in humans, Craig [2009, 2010] proposed 
that this cortical territory in the human dorsal anterior 
insula in may lack an obvious morphological homologue in 
the monkey. On the basis of the extensive evidence from 
the functional neuroimaging literature, it was shown that 
the human dorsal anterior insula has a crucial role in inte- 
grating salient information from all neural systems at the 
present instant, which suggests that in humans, "right and 
left anterior insulae are crucial components of a network 
that engenders human awareness" [Craig, 2010]. Although 
it is difficult on the basis of the current evidence to deter- 
mine whether the anterior insula of humans represent a 
newly evolved structure, we mention that one allometric 
study [Semendeferi and Damasio, 2000] reported that the 
human insula may have increased more than expected dur- 
ing hominid evolution; however, the results did not reach 
significance, and the authors concluded that the amount of 
examined MR images (N = 29, 10 humans) was not large 
enough to be confident in the obtained results. 

Suggested Asymmetries in the 
Tractography Results 

Our tractography results suggested differences in the 
connectivity maps of the right (RH) and left (LH) insula. 



We therefore performed additional analyses to study inter- 
hemispheric differences in the tractography samples to the 
target regions chosen to perform connectivity-based classi- 
fication, and we considered differences in the extent of the 
insula connected to each target, and the variability across 
subjects of the centers of gravity of the connected insular 
territory. To our knowledge, evidence of lateralization of 
insular connections in macaques is lacking; however, 
functional neuroimaging investigations in humans added 
valuable information about functional and structural asym- 
metries in the insula and provide in some cases arguments 
for a discussion of our tractography findings. Watkins 
et al. [2001] analyzed 142 healthy adults with voxel-based 
morphometry and found RH > LH asymmetry in the 
amount of gray matter in the anterior inferior insula. In a 
study on 300 MRI images and 16 fixated brain specimen, 
Naidich et al. [2004] documented the variability in the con- 
figuration of insular gyri as well as differences in this con- 
figuration between LH and RH insulae, often in the same 
subject. A recent meta-analysis on 1768 fMRI experiments 
suggests the presence of asymmetries in the insula for 
some functional categories such as taste, interoception, 
somatosensation, empathy, and speech [Kurth et al., 
2010b]. On the basis of a great amount of evidence indicat- 
ing a forebrain emotional asymmetry, Craig [2005, 2009, 
2010] proposed that RH and LH anterior and mid insula 
are associated with sympathetic and parasympathetic activ- 
ity, respectively, and that their coordinated activity 
embodies a neural system aimed at optimizing energy ex- 
penditure. Finally, we mention that a recent work used 
metrics from graph theory on whole-brain tractography 
and found a RH > LH asymmetry for measures of global 
efficiency and interconnectivity, as well as for regional 
measures of betweenness centrality, or relative importance 
of several brain regions for the connectivity of the whole 
network [Iturria-Medina et al., 2011]. Results in the present 
work suggest on a global level that the RH anterior insula 
has a more widespread connectivity patterns with respect 
to the LH anterior insula (see Fig. 4). Although the results 
from Iturria-Medina and colleagues did not find significant 
difference of betweenness centrality for the RH or LH insu- 
lae, we suggest the possibility that such difference could be 
found by comparing only the anterior insulae. 

The insular region sending above-threshold samples to 
the Limbic ROI was found to be significantly bigger in the 
LH, even though the difference in the number of samples 
recorder in the Limbic ROI did not significantly differ 
between hemispheres. It is possible that the detection of 
connectivity with the limbic lobe was favoured in the LH 
by a reported leftward asymmetry of fractional anisotropy 
(LH > RH) of the subinsular portion of the uncinate fasci- 
culus [Kubicki et al., 2002; Rodrigo et al., 2007]; on the 
other side, a completely different hypothesis points to the 
presence of a region in the mid-posterior ventral insula, 
consistent with the location our results (see Fig. 7) which 
was found in a recent meta-analysis to be specifically 
involved in empathy [Kurth et al., 2010b]. 
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Samples to the primary somatosensory cortex (SI) were 
sent from the middle insula in both hemisphere; however, 
the RH insula was characterized by a second territory of 
consistent projections across subjects, located in the poste- 
rior insula. An interesting hypothesis is that this asymme- 
try is related to the extensively described RH posterior 
insula involvement in the sense of limb ownership in both 
healthy subjects and patients (see Ibanez et al. [2010] and 
Karnath and Baier [2010] for reviews). Several clinical 
studies [Baier and Karnath, 2008; Cereda et al., 2002; Kar- 
nath et al., 2005] reported that in patients with hemipare- 
sis/hemiplegia, the presence of anosognosia and other 
"disturbed sensation of limb ownership" (DSO [Baier and 
Karnath, 2008]) is more frequent in cases of damages to 
the RH posterior insula. In healthy subjects, the rubber 
hand illusion paradigm was used to manipulate the sense 
of limb ownership by Tsakiris et al. [2007]. In this PET 
experiment, participants saw either a right or a left hand 
being rubbed, in or out of synchrony with their own hid- 
den right hand; they had to indicate the position of their 
limb before and after the stimulation. A correlation was 
observed between brain activity and the proprioceptive 
drift toward the rubber hand, in a region of the posterior 
RH insula that co-localizes with the area of maximum 
overlap in the posterior RH insular cluster found to be 
connected to SI (see Fig. 4A in Tsakiris et al. [2007]). 

Several insights regarding the hemispheric asymmetry 
of insular connectivity with BA44 come from the extensive 
functional neuroimaging literature on the lateralized 
involvement of the insula in tasks related to speech and 
singing. For example, in the speech domain, the RH 
fronto-opercular region is involved in extracting slow pitch 
information from the speech stream, while the LH fronto- 
opercular region performs extraction of segmental units 
from the speech signal [Meyer et al., 2002; Wildgruber 
et al., 2004]. This evidence is strikingly in line with the ob- 
servation that LH insular activity increases and the RH in- 
sular activity decreases, as a function of the rate of stimuli 
presentation (trains of clicks and syllables) [Ackermann 
and Riecker, 2004; Ackermann et al., 2001]. A lateralization 
effect is also suggested from experiments on the different 
meaning conveyed by prosody during speech production: 
while production of linguistic prosody (e.g., formulating a 
question) is associated with activity in the bilateral IFG 
and the LH insula, emotional prosody (e.g., a happy into- 
nation) is associated with activity in the left IFG and in the 
RH insula [Aziz-Zadeh et al., 2010]. Finally, in a task 
involving either overt recitation of the months of the year 
or overt singing a nonlyrical tune, Riecker reported the 
same asymmetry for both insula and IFG: LH activity in 
these regions was associated with the speech production 
task, while RH activity was associated with the singing 
task [Riecker et al., 2000]. Given these evidence, and con- 
sidering in general the role of the insula in dynamically 
integrating salient aspects of inner and outer environment 
to generate appropriate energy-efficient behavioral 
responses [Craig, 2009; Mesulam and Mufson, 1985], we 



suggest that the stronger connectivity between BA44 and 
RH insula found in our tractography study would make 
the right hemisphere better suited to convey the emotional 
state of the organism during speech production [Acker- 
mann and Riecker, 2004] and, in general, during social 
interaction as well as during singing. Also, the RH insula 
could be particularly well suited for the on-line integration 
of somatosensory, kinesthetic, and auditory feedbacks, 
which is required for the enhanced vocal motor control 
during singing [Kleber et al., 2010]. 



Intersubjective Variability 

In the present study, we also considered the intersubjec- 
tive variability in the connection probability of the insula 
with several brain regions. We reported, for each target 
mask and hemisphere, spatial variability maps of the COG 
of the connected insular cluster across subjects (Fig. 7, up), 
quantified the spatial spread as a function of hemisphere 
and target, and computed consistency maps showing for 
each target mask how many subjects sent an above-thresh- 
old number of samples from each insular voxel (Fig. 7, 
up). This analysis suggested a considerable spatial vari- 
ability, particularly in the right hemisphere, both in the 
COG and in the location of the insular voxels found to be 
connected with BA6, SI, and the Limbic ROI. Regarding 
the connectivity between insula and Limbic ROI, which 
was observed especially in the ventralmost anterior insula, 
we point out that this region is characterized by the pres- 
ence of the transverse insular gyrus (TG). This region 
extends from the insular pole and apex to the posterome- 
dial orbital lobule [Naidich et al., 2004; Ture et al, 1999] 
(see Fig. 2 in Naidich et al. [2004]) and represents the loca- 
tion where the densest concentration of VEN in the insula 
was found in humans (see Allman et al. [2010]; Fig. 3). 
Because this region is phylogenetically young [Rose, 1928], 
Craig [2010] proposed that large interindividual structural 
differences might be observed in the TG, and various dis- 
section studies indeed reported considerable morphologi- 
cal variability in this region [Naidich et al., 2004; Ture 
et al., 1999; Varnavas and Grand, 1999]. 

Various studies on postmortem specimens and MR 
images reported a considerable variability in the sulcal 
and gyral morphology of the insular cortex in humans 
[Afif and Mertens, 2010; Johnson et al, 2009; Naidich 
et al., 2004; Ture et al., 1999; Varnavas and Grand, 1999] as 
well as across species [Butti and Hof, 2010]. The intersub- 
jective variability in the spatial distribution of different 
cortical types was reported in a recent cytoarchitectonic 
parcellation study on the postcentral insular region [Kurth 
et al., 2010a]; however, in that study, the relationship 
between cytoarchitecture and sulcal variability was not 
examined. Intersubjective variability in sulcal morphology 
represents a well-known problem for group analyses in 
neuroimaging, because a misregistration of the sulci across 
subjects contributes to increase the probability of 
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performing group analyses on nonhomologous brain 
regions. Fully aware of these issues, when comparing the 
relative rate of change of the connectivity patterns in the 
insula and in the mPMC, we calculated the relevant meas- 
urements in each subject, hemisphere, and region of inter- 
est separately. The metric of interest was represented by 
the magnitude of the gaps in each subject's LE, with larger 
gaps denoting a more abrupt transition between connectiv- 
ity patterns. By comparing the LE of insula and mPMC of 
each subject, Wilcoxson signed rank test reported signifi- 
cantly larger gaps in the LE of the mPMC, suggesting that 
the transition between different connectivity patterns is 
more gradual in the insula than in the mPMC. Because the 
gaps are calculated in the LE of each subject, hemisphere, 
and region of interest separately, rather than on a group 
average, the recovered difference in the rate of change of 
connectivity patterns cannot be affected by issues due to 
averaging across subjects, such as smoothing or misregis- 
tration across different brains. Regarding the issue of sul- 
cal variability, we note that this analysis compares the 
magnitude of the gaps in the space of the LE of each sin- 
gle subject, independently of the anatomical location of 
those gaps, and therefore is not affected by intersubjective 
differences in sulcal morphology. 

In this work, we report maps of sulcal variability in our 
sample (see Fig. 7, down), although here we do not inves- 
tigate the link between the differences in the localization 
of the trajectories of connectivity variability and macroana- 
tomical landmarks such as sulcal patterns. In which way, 
the morphological variability of the insula relates to 
cytoarchitecture, and functional difference is an important 
topic for future multimodal imaging studies. We note 
however that the available literature about humans do not 
reflect a general consensus regarding the relationship 
between cytoarchitecture and the sulcal anatomy (see also 
Supporting Information Fig. S5). In 1909, Brodmann di- 
vided the insula in two main cortical fields and stressed 
that their border, although in continuation with the central 
(Rolandic) sulcus, did not correspond exactly with the cen- 
tral insular sulcus [Brodmann, 1909; 122). The parcellation 
presented by Rose [1928] also displays some relationships 
with the sulcal morphology, although in some cases, and 
specifically in the middle and posterior short gyri, differ- 
ent cortical types are present. Even more complex is the 
scheme presented by Brockhaus [1940], in which the dif- 
ferentiation between the main cortical types is almost or- 
thogonal with respect to the sulcal morphology. Similarly, 
the study performed by Bonthius et al. [2005] is in line 
with the report of Mesulam and Mufson [1982a] about the 
radial organization of the cytoarchitecture. Finally, the 
recent parcellation study by Kurth et al. [2010a] on the 
postcentral insular cortex maps the granular territories of 
the insula (Igl and Ig2) in the dorsal part of the postcen- 
tral gyri, caudally located with respect to the central sul- 
cus as well as a dysgranular sector (Idl) located mainly in 
the posterior part of the IPS. The latter study is particu- 
larly important in the context of clarifying the relationship 



between sulci and cytoarchitecture, because it is the first 
observer-independent parcellation study of the insula. 
Maps of the anterior sector of the insula obtained with this 
methodology will help clarify whether a stronger consis- 
tency between cortical types and sulcal morphology can 
be recovered in the precentral insular territories. 



A Comment on the Claustrum 

The gray matter of the claustrum is located underneath 
the insular cortex. Middle and dorsal insular territories are 
separated from the claustrum by the extreme capsule, 
while the inner stratum of the agranular insula is coexten- 
sive with the claustrum [Bonthius et al., 2005; Mesulam 
and Mufson, 1982a]. Given the close proximity of insula 
and claustrum, and the spatial resolution of the acquired 
DW images, it is possible that some claustral gray matter 
was included in the ventral part of the insular ROI used to 
perform probabilistic tractography. One effect of this inclu- 
sion could have been the presence of samples located in 
extrastriate regions of the occipital lobe (lingual gyrus and 
lateral putative BA19), because, to our knowledge, no 
direct connections between the insula and the occipital 
lobe have been reported in animal studies. On the other 
hand, animal studies have shown the direct involvement 
of the claustrum in processing visual stimuli (see referen- 
ces in Crick and Koch [2005], Mufson and Mesulam 
[1982], Tanne-Gariepy et al. [2002]) as well as claustral 
connectivity with parieto-occipital regions [LeVay and 
Sherk, 1981; Shipp et al., 1998]. Moreover, a recent com- 
bined dissection and DTI study on the claustrum in 
humans [Fernandez-Miranda et al., 2008] suggested the ex- 
istence of white-matter pathways between the claustrum 
and the occipital lobe via the inferior fronto-occipital fasci- 
culus proposed by Catani et al. [2002] (but see a discussion 
on the existence of this fiber bundle in Schmahmann and 
Pandya [2007] and Schmahmann et al. [2007]). 

If this partial volume effect would have had a major 
influence on the recovered structure of connectivity vari- 
ability, our results should have displayed some consis- 
tency with the available knowledge from animal studies 
about the patterns of histological, functional, and /or con- 
nectivity variation in the claustrum that are known to dif- 
fer from those in the insula [Brand, 1981; Kowianski et al., 
1999; Mufson and Mesulam, 1982]. In contrast to the 
insula, whose connectivity is organized mainly along a 
rostrocaudal gradient [Mesulam and Mufson, 1982b] with 
its various territories showing comparatively distinct con- 
nectivity patterns, the claustrum is divided histologically 
[Brand, 1981] and morphometrically [Kowianski et al., 
1999] in a dorsal and a ventral section, with comparatively 
unspecific projection to different cortical regions [Crick 
and Koch, 2005; Tanne-Gariepy et al., 2002]. It is therefore 
reasonable to conclude that the effect of potentially includ- 
ing part of the claustrum in our seed region must have 
had at most a minor effect on our results. 
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Gradual Variation of Connectivity Patterns 

Connectivity-based parcellation has been previously 
used to identify clusters with distinct connectivity in the 
medial premotor wall [Johansen-Berg et al., 2004; Klein 
et al, 2007], in the IFG [Anwander et al., 2007; Klein et al., 
2007], in the cingulate cortex [Beckmann et al., 2009], and 
in the lateral premotor cortex [Tomassini et al., 2007], 
whose boundaries were consistent with available cytoarch- 
itectonic atlases and with results from animal studies. In 
these studies, spectral reordering [Barnard et al., 1995] or 
k-means clustering [Hartigan and Wong, 1979; Hartigan, 
1975] were used, because a sharp segregation between 
clusters with distinct connectivity was expected on the ba- 
sis of previous evidence from tracer-injection studies in 
macaques [Behrens and Johansen-Berg, 2005]. Accordingly, 
the relevant question was where the boundaries separating 
different clusters are located in humans. However, hard 
clustering cannot examine whether it is appropriate to 
assume that the data contain distinct clusters. In the case 
of the insula, several animal studies suggest a gradual 
rather than abrupt variation in connectivity: the gradient 
of increasingly more complex architecture spanning from 
agranular to granular territories [Friedman et al., 1986; 
Jones and Burton, 1976; Mesulam and Mufson, 1982a; 
Mesulam and Mufson, 1985; von Economo, 1927] appears 
to be related to the gradual change in connectivity pat- 
terns, as reported by tracer-injection studies [Amaral and 
Price, 1984; Chikama et al., 1997; Friedman et al, 1986; 
Fudge et al., 2005]. On other hand, a very recent study on 
the cytoarchitecture of the posterior human insula 
reported the presence of clear boundaries between granu- 
lar and dysgranular insular regions [Kurth et al., 2010a]. 
Neuroimaging studies are only beginning to be able to 
detect differences among different cytoarchitectonic fields 
in vivo [Eickhoff et al, 2005a; Fatterpekar et al., 2002; 
Geyer et al., 2011], and, in the present study, we were not 
in the position to investigate the relative smoothness of 
the transition across cortical types in the insula. We 
instead used the results from probabilistic tractography (1) 
to investigate the likelihood of either a sharp (clustered) or 
a gradual transition among connectivity patterns and (2) 
to compare the relative sharpness of this transition with 
another brain region (mPMC) where sharply segregated 
clusters of connectivity patterns were shown in previous 
studies. 

In a previous study using repeated fc-means clustering, 
we showed that applying hard clustering to the cross-cor- 
relation matrix of connectivity patterns is undesirable in 
the insula: repeated application of fc-means clustering 
reveals that in comparison with the mPMC, a significantly 
higher proportion of the middle insula cannot be reliably 
categorized using this method [Nanetti et al., 2009]. In this 
work, we therefore applied LE [Belkin and Niyogi, 2002] 
in order to analyze how the connectivity of the insula 
varies across different territories. Given the connectivity 
map of each seed voxel, which is the output of probabilis- 



tic tractography, this algorithm uses the information about 
the similarity between each and every other connectivity 
map to recover, if present, the main trajectories of connec- 
tivity variation across different territories of the seed 
region. In the case of the insula, several anatomical studies 
suggested that the main trajectory of connectivity variation 
would be topographically localized along the rostrocaudal 
anatomical axis and that the variation across connectivity 
patterns of different insular territories would be gradual 
rather than abrupt. We therefore expected LE to represent 
the change in connectivity patterns as a one-dimensional 
structure with relatively small interruptions, especially 
compared to other cortical regions where the change in 
connectivity is more sharp. In the latter case, we expected 
to detect relatively bigger gaps in the structure of connec- 
tivity variability recovered by LE. Note that while prevent- 
ing artifactual clustering not supported by the data, this 
algorithm is nonetheless capable of revealing clusters 
when present [Belkin and Niyogi, 2002] and does not 
require an a priori hypothesis about the number of clus- 
ters. Additional precautions were also taken to avoid that 
the results of LE could be affected by the mere anatomical 
distance within seed voxels as well as between seed and 
target voxels. Specifically, (1) we corrected the count of the 
samples in each target voxel by the mean length of the 
pathways connecting each seed with each target voxel to 
ensure that the results of tractography would have not 
been biased toward anatomically close target voxels. The 
same correction was used for this reason in connectivity- 
based parcellation/classification studies [Beckmann et al., 
2009; Eickhoff et al., 2010; Tomassini et al., 2007]. (2) An 
Euclidean anatomical distance matrix was added to the 
correlation/similarity matrix of the connectivity maps to 
penalize the similarity between seed voxels, which are 
very close in the anatomy. This prevents that the connec- 
tivity maps of two seed voxels would be quantified as 
similar in the LE just because of their anatomical proxim- 
ity. (3) Also, potentially outlier points produced by an ini- 
tial LE were removed prior to any analysis of the rate of 
change of connectivity to avoid the possibility that conti- 
nuity could be suggested by the presence of noise (see 
Supporting Information: Theory and implementation of 
LE). 

Applying this methods to the results of probabilistic 
tractography from the insula revealed a trajectory of vari- 
ability between two extremes featuring very different 
connectivity maps. Examining the spatial location of this 
trajectory on the insular surface revealed that these 
extremes were located in the anterior dorsal and ventral 
insula and in the dorsal and posterior part of the long in- 
sular gyri, respectively. Visual inspection of the LE sug- 
gested that this change of connectivity was gradual 
rather than abrupt. To assess this graduality quantita- 
tively, we compared this profile with that obtained by 
applying the same technique to the mPMC. We found 
that LE also revealed a one-dimensional structure of con- 
nectivity change in the mPMC, corresponding the the 
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rostrocaudal anatomical axis. However, the change of 
connectivity was significantly more gradual in the case of 
the insula, while that in the mPMC was characterized by 
the presence of larger "gaps" that are consistent with the 
results of hard clustering methods (see Fig. 5) [Anwander 
et al., 2007; Johansen-Berg et al., 2004; Klein et al., 2007; 
Nanetti et al., 2009]. These results contribute to our 
understanding of the connectivity of the insula in 
humans and suggest that LE (and possibly other methods 
for nonlinear dimensionality reduction, see Saul et al. 
[2005] for a review) can be useful in analyzing the con- 
nectivity structure of different regions in the brain, 
including those where the transition between different 
types of connectivity is not sharp. 

Given the present results as well as the evidence from 
animal studies, it comes natural to ask how a gradual tran- 
sition in connectivity could relate to the integrative func- 
tionality of the insular cortex. In a classical work, 
Mesulam and Mufson noted that "the insula as well as the 
adjacent orbitofrontal and temporopolar areas are charac- 
terized by a remarkable heterogeneity in cortical architec- 
ture, connectivity, and physiology. This heterogeneity 
which is a hallmark of paralimbic (mesocortical) parts of 
the brain, indicates neither a chaotic organization nor a 
haphazard distribution of behavioral specialization. 
Instead, evidence indicates that the insula, like other com- 
ponents of the paralimbic brain, acts as a multipurpose 
region of cortex were a remarkably wide range of neural 
processes modulate behaviors which primarily depend on 
interactions between the extrapersonal world and the in- 
ternal milieu" [Mesulam and Mufson, 1985]. These consid- 
erations were substantialized during the following decades 
of functional neuroimaging works: recently, Craig [2009], 
while detailing his model of human awareness, proposed 
that a complete representation of all salient information 
for the organism at a given time is produced in a poste- 
rior-to-anterior progression in the insula by means of a 
progressive integration of "homeostatic conditions with 
the sensory environment, and with motivational, hedonic 
and social conditions in other parts of the brain" [Craig, 
2009]. In a recent meta-analysis, Kurth et al. [2010b] 
focused on both functional integration and functional spe- 
cialization in the insula: in all the examined functional do- 
main, the regions that showed a specific effect to a certain 
functional category were substantially smaller than the 
regions involved in multiple different categories. The 
authors comment this evidence as an indication that the 
insula is characterized both by functional specialization 
and functional integration. Interestingly, several of the 
regions identified with the meta-analysis of functional inte- 
gration span different territories featuring different 
cytoarchitecture, according to the available maps, and, 
according to our results, progressively different connec- 
tional fingerprints. In the same work, Kurth and col- 
leagues propose that the presence of intrainsular 
connections and the evidence of overlap between func- 
tional categories in the anterior dorsal insula could be two 



possible ways of reasoning about how functional integra- 
tion could be embodied in the insula. It is outside the 
scope of this work to define the most plausible or efficient 
hodological model for functional integration in the cortex; 
however, we speculate that the relatively gradual shift in 
connectivity patterns, which was recovered in the insula, 
would be consistent with the concept of a cortical region 
whose nature is fundamentally integrative, because such a 
region should have the capability of both receiving differ- 
ent sources of information (interoceptive, sensory, emo- 
tional, and cognitive), and the different aspects of them 
(for instance the physical and affective components of 
touch), abstracting them from their original domain and 
finally integrating them according to their saliency for the 
present situation that the organism has to cope with. Sharp 
functional subdivisions are usually accompanied by sharp 
transition of connectivity, like in the mPMC, while the 
integrative role of a region such as the insula could be bet- 
ter substantialized by the close proximity of different terri- 
tories with slightly different, nevertheless similar or 
related, connectivity and functionality, sensitive to differ- 
ent aspects of the processed information. This is indeed 
the situation suggested on a regional scale by the present 
findings obtained with LE, at the spatial resolution 
allowed by in vivo DWI, as well as by previous studies in 
macaques, and in particular regarding insular connections 
with the amygdala [Amaral and Price, 1984] and the basal 
ganglia [Chikama et al, 1997; Fudge et al., 2005]. The 
same situation is suggested when looking at the segrega- 
tion of whole brain projections from either the mPMC or 
the insula: SMA and preSMA (F3 and F6) contain mark- 
edly different connectivity with the rest of the brain and 
relatively modest connections among each other [Luppino, 
1993]; on the other side, as we highlighted in the introduc- 
tion, dysgranular territories of the insula project to regions 
in the brain that are also targeted by agranular and granu- 
lar territories, and their differential connectivity is reflected 
in their topographical proximity with either agranular or 
granular insular territories. 

Comparison with Functional Studies 

Insula and speech articulation 

Many studies reported an involvement of the left ante- 
rior insula in overt language production tasks and specifi- 
cally in the motor control of speech [Ackermann and 
Riecker, 2004, 2010; Bamiou et al, 2003; Blank et al., 2002; 
Duffau et al., 2000; Price, 2000; Riecker et al., 2005; Saur 
et al., 2006; Wise et al., 1999]. Deficits in speech coordina- 
tion are associated with lesions [Dronkers, 1996] and hypo- 
metabolism [Nestor et al., 2003] in the short gyri of the 
insula, and its activity increases as a function of syllabic 
complexity and sequence complexity during a polysyllabic 
items production task [Bohland and Guenther, 2006]. In 
the perceptual domain, the activity of the left anterior 
insula is positively (left hemisphere) or negatively (right 
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Comparison with functional studies. The average trajectory of 
connectivity variability is overlaid onto the MNI single subject 
template. Only voxels belonging to at least half of the partici- 
pants (five subjects) are mapped. Blue dots indicate the y and z 
coordinates of the local maxima for fMRI studies on motor con- 
trol of speech in the following studies: (I) Dronkers [1996], (2) 
Blank et al. [2002], (3) Nestor et al. [2003], (4) Plante et al. 
[2006], (5a) and (5b) Saur et al. [2006], and (6) Geiser et al. 
[2008]. These locations are all in the insular territories where 
the strongest projections to BA44, BA45, and BA6 were found 
(see Fig. 6). Red and yellow dots indicate the locations from 
Ostrowsky et al. [2002] where intracortical stimulations elicited 
either painful (yellow) or nonpainful (red) sensations. Note the 



prevalent location in the long gyri of the insula, which we found 
to have maximum connection probability with somatosensory 
and posterior parietal regions (see Fig. 6). The meta-analysis of 
Schweinhardt et al. [2006] identified in the rostral anterior 
insula (rAI), caudally delimited by the short insular sulcus 
(dashed black line), as specifically involved in clinical pain with 
respect to nociceptive (experimental) pain, involving more pos- 
terior insular regions. The limbic lobe was found in the present 
study to receive projections mostly from the same insular terri- 
tories (anterior short gyrus and ventral anterior insula — see Fig. 
6). [Color figure can be viewed in the online issue, which is 
available at wileyonlinelibrary.com.] 



hemisphere) correlated with the presentation rate of trains 
of clicks or syllables [Ackermann et al., 2001], while this 
effect is not present in patients with developmental dys- 
lexia [Steinbrink et al., 2009]. Also, both the left anterior 
insula and BA44/45 in the IFG are involved in detecting 
pure rhythmic aspects of speech production [Geiser et al., 
2008] as well as in recognizing degraded statements on the 
basis of their prosodic profile [Plante et al., 2006]. The fact 
that these anterior regions of the insula are more involved 
in language functions is compatible with the fact that this 
insular territory had the highest connection probability 
with the key frontal language regions (BA44, BA45, and 
BA6; see Fig. 6) in our study. 

Insula and pain 

The insula has been considered to be a core region for 
the integration of affective, empathic, interoceptive, and 
somatosensory components of pain [Brooks and Tracey, 
2007; Lamm and Singer, 2010] with the insula being the 
most frequently reported brain region in pain-related fMRI 
experiments (see Apkarian et al. [2005] for a review). Two 
themes emerge from this vast number of publications: first, 
the rostral short gyri of the insula (rAI) seem to be most 
involved in the affective component of pain [Hein and 
Singer, 2008; Schweinhardt et al, 2006; Singer et al., 2004]; 
second, more posterior territories of the insula represent 
the sensory aspects of pain, with activity in these regions 



correlating with the physical properties of the stimulus 
[Frot et al., 2007] and electrostimulations determining sen- 
sory (tactile or nociceptive) qualia [Ostrowsky et al., 2002] 
(see Fig. 8). These findings are compatible with our results. 
First, in our study, we recorded the strongest connectivity 
to the limbic ROI (amygdala, EC, and hippocampus) in 
the most ventral anterior insula, and, to a lesser extent, in 
the entire anterior short gyrus and in the middle ventral 
insula. Second, the regions we found to be most connected 
with somatosensory and posterior parietal cortices were 
located in the posterior short gyrus and in the long insular 
gyri, which correspond to the regions found to be mostly 
involved in processing the physical properties of the stim- 
ulus (see Fig. 8). 

Connection probability and functional 
organization of the insular cortex 

In addition to the above considerations regarding speech 
articulation and pain, we report the results of two recent 
meta-analyses [Kurth et al., 2010b; Mutschler et al, 2009] 
performed on results from functional neuroimaging stud- 
ies, in order to show the relationship of those findings 
with the results of our connectivity-based classification. 

Mutschler and colleagues [Mutschler et al., 2009] per- 
formed a meta-analysis of the functional organization of 
the anterior insula using the results of 58 studies reporting 
149 coordinates. The considered studies included auditory 
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and language processing, motor tasks involving upper and 
lower limbs, responses related to peripheral autonomic 
changes as well as tasks resulting in co-activation of the 
insula with the amygdala. The resulting maps are shown 
in Figure 6 (B1-B3) to appreciate the topographical rela- 
tionship between functional specialization and peaks of 
connection probability. The functional meta-analysis map 
considering peripheral and autonomic changes as well as 
co-activation with the amygdala is localized in the anterior 
ventral insula, very close to the region where we found 
the peak for the connection probability with the limbic 
ROI (encompassing amygdala, EC, and hippocampus). The 
functional map considering studies on auditory and lan- 
guage processing is completely enclosed in the insular 
regions, where maximum connection probability was 
found with BA44 and BA45. Finally, the peak of the meta- 
analysis for studies on motor activity was located in a spe- 
cific territory in the middle insula, in the ventral part of 
the posterior short gyrus, and to a lesser extent in the mid- 
dle short gyrus, in the same locations, where we found 
connection probability values with the premotor cortex 
(BA6). 

Kurth et al. [2010b] performed a meta-analysis on the 
coordinates obtained from 1768 experiments and collected 
in the BrainMap database [Fox et al., 2005b; Laird et al., 
2005, 2009] and Pubmed. In this study, the authors investi- 
gated both functional integration and functional differentia- 
tions in the insula in 13 different functional categories (see 
Kurth et al. [2010b], Table 1). When combining these catego- 
ries into four functional domains (sensorimotor, cognitive, 
chemical sensory, and social-emotional), the insula dis- 
played a domain-specific functional organization, which 
mostly matches our results obtained by connectivity-based 
parcellation, in particular when relating BA44 and BA45 to 
the language-cognitive domain, the connection probability 
with the Limbic ROI to the social-emotional domain, and 
the connection probability with BA6 and SI to the somato- 
sensory domain. In sum, these comparisons increases our 
confidence in the results obtained by performing in vivo 
probabilistic tractography of the insula and supports the hy- 
pothesis that the differential connectivity and the differen- 
tial cytoarchitecture of the insula represents the substrate of 
the functional specialization of different insular territories. 

Parcellation of the insular cortex on the basis 
of resting state functional connectivity 

Recently, two works [Cauda et al., 2011; Deen et al., in 
press] used resting resting-state functional MRI (RS fmri) 
to perform a voxelwise parcellation of the insula on the 
basis of the correlation between the time course of each in- 
sular voxel with that of each other voxel in the brain. This 
methodology aims at recovering the intrinsic functional 
architecture of the brain by looking at regions, which fea- 
ture similar spontaneous fluctuations of the blood-oxigen- 
ated level dependent signal. In the study of Deen et al. [in 
press], fc-means clustering was performed on the average 



correlation matrix of insular functional connectivity pro- 
files across subjects, and the optimal number of k was set 
to 3. Two clusters were identified in the anterior insula, 
the third in the posterior insula. In the study of Cauda 
et al. [2011], fuzzy c-means clustering was performed for 
each subject independently, and the optimal number of 
clusters was quantitatively set to 2. The algorithm divided 
the insula in one anterior and one posterior cluster. 
Although the spatial consistency of the anterior cluster 
found in Cauda et al. [2011] with the two anterior clusters 
found by Deen et al. [2010] cannot be precisely assessed 
on the basis of the published data as well as the location 
of the recovered subdivision with respect to the sulcal 
anatomy, we note that in the study of Cauda et al. [2011], 
the region where an overlap between the two clusters was 
found is located approximately where our previous 
repeated fc-means clustering study on tractography data. 
Nanetti et al. [2009] found the highest density of voxels 
that could not be reliably assigned to either clusters: in the 
middle and posterior short insular gyri and in the anterior 
part of the long insular gyri. This convergent result from 
two different imaging modalities suggests that voxels in 
this insular region feature a high uncertainty when cluster- 
ing techniques are used for the purpose of parcellation. 

The discrepancy in the optimal number of clusters 
found between the study of Cauda et al. [2011] and that of 
Deen et al. [2010] might be due to a number of factors. For 
instance, Cauda et al. [2011] used spatial smoothing to to 
increase local correlation across neighbouring voxels, 
which could have masked the differences between dorsal 
and ventral anterior insula, as found by Deen et al. 
[2010]; on the other side, Deen et al. [2010] averaged all 
the individual subjects' correlation matrices before per- 
forming fc-means clustering, while performing clustering 
on the single-subject correlation matrices would have 
made possible to check the consistency of the spatial loca- 
tion of clusters in order to decide the appropriate number 
of k (as it was done by Beckmann et al. [2009]). Further 
investigations by means of other dimensionality reduction 
techniques, which do not require an a priori choice of the 
number of clusters, such as multidimensional scaling or 
also the LE used in this study, could also help to clarify 
this point. In addition, these techniques could give further 
insights regarding whether the region of uncertainty 
recovered by Cauda et al. [2011] reflects a situation where 
the transition between functionally different cortical terri- 
tories is less abrupt than in other brain regions [Kim et al., 
2010]. 



CONCLUSIONS 

In vivo probabilistic tractography is a relatively recent 
field of research, which still suffers from many limitations 
including low-signal-to-noise ratio in the images, low spa- 
tial resolution, problems related to tracking fibres in 
regions of particular structural complexity. In addition, it 
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cannot discriminate between afferent and efferent projec- 
tions, and the recorded number of samples in the target 
voxels can be influenced by the distance between seed and 
target as well as by the geometry of the pathways connect- 
ing seed and targets [Johansen-Berg and Behrens, 2006; 
Jones, 2008]. Because of these limitations, the golden 
standard for investigations in anatomical connectivity is 
still represented by tracer-injection studies in animals, and 
it is therefore wise at the moment to focus current proba- 
bilistic tractography studies on testing in humans specific 
hypotheses derived from animal evidence. Considering the 
macaque data as the basis for generating hypothesis and 
as a rough meter of the validity of the findings reduces 
the likelihood of errors. Where the two sources of informa- 
tion diverge, a species difference or a methodological limi- 
tation might exist. Disentangling those two possibilities is 
difficult at the current stage of methodological develop- 
ment. Keeping these caveats in mind, we used probabilis- 
tic tractography to examine for the first time the 
connectivity of the human insula in vivo. The hypotheses 
that guided our work were grounded on the evidence 
about a relation between connectivity and cytoarchitecture 
in macaques and about the similar topographic arrange- 
ment of the cytoarchitecture in human and macaque 
insula. Fiber tracking was performed with an algorithm ca- 
pable of resolving two main fiber populations within a 
voxel [Behrens et al., 2007], and the number of tractogra- 
phy samples reaching each target voxel was corrected by 
the mean path length between seed and target, as in previ- 
ous studies [Beckmann et al., 2009; Eickhoff et al., 2010; 
Tomassini et al., 2007]. By analyzing the results with LE, 
we recovered in the insula a trajectory of connectivity vari- 
ability spanning the rostrocaudal anatomical axis. Also, 
the transition among connectivity patterns in the insula 
was significantly smoother than that in the mPMC. The 
brain regions reached by tractography from different insu- 
lar territories were largely consistent with the available 
evidence from animal studies both in terms of specificity 
(i.e., lack of aberrant connections) and sensitivity (i.e., find- 
ing those connections known to exist), and further analysis 
of the tractography samples count in target voxels, as well 
as of the connection probability with cytoarchitectonically 
defined regions suggested the presence of a relationship 
between connectivity patterns and cytoarchitecture, as it 
was found to be the case in macaques. Preliminary results 
from our laboratory suggest that very similar results can 
be obtained by performing the same analyses on ex vivo 
high-resolution diffusion-weighted data from macaque 
(see Cerliani [2009] and Supporting Information Fig. S2). 

Although the available in vivo tractography techniques 
are still far from achieving results comparable to those 
derived from invasive studies, our study represents a fur- 
ther confirmation that tractography techniques can be 
used to test specific anatomical hypotheses — such as the 
relative sharpness /smoothness of a transition between dif- 
ferent connectivity patterns — when these can be formu- 
lated at a degree of resolution compatible with the current 



limitations of DW imaging and white-matter tractography. 
More generally, we hope that our results will significantly 
broaden the spectrum of tools to analyze the results of in 
vivo tractography. The hard-clustering methods that have 
been used so far required the experimenter to think about 
the anatomical configuration of a region in terms of dis- 
tinct clusters with sharp borders. This limited the applica- 
tion of in vivo connectivity-based parcellation studies to 
structures where such clusters are to be expected. LE, on 
the other hand, provides the opportunity to first examine 
whether a region's connectivity pattern is best described in 
terms of clusters featuring sharper transitions, as for the 
mPMC, or in terms of a more gradual transition, as for the 
insula. In the former case, Laplacian eigenmaps are useful 
to help estimating the number of clusters to ask to a clus- 
tering algorithm such as fc-means. In the latter, it affords 
the possibility to describe the differences and similarities 
among the connectivity patterns of various territories of a 
brain region, without artificially dividing it into sharply 
segregated subregions, which could provide a suboptimal 
or incomplete representation of the intrinsic organization 
of the data. 
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